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ABSTRACT 

We have observed the Virgo Cluster spiral galaxy, NGC 4845, at 1.6 and 6 GHz using the 
Karl G. Jansky Very Large Array, as part of the ‘Continuum Halos in Nearby Galaxies - an 
EVLA Survey’ (CHANG-ES). The source consists of a bright unresolved core with a surrounding 
weak central disk (1.8 kpc diameter). The core is variable over the 6 month time scale of the 
GHANG-ES data and has increased by a factor of « 6 since 1995. The wide bandwidths of 
GHANG-ES have allowed us to determine the spectral evolution of this core which peaks between 
1.6 and 6 GHz (it is a GigaHertz-peaked spectrum source). We show that the spectral turnover 
is dominated by synchrotron self-absorption and that the spectral evolution can be explained by 
adiabatic expansion (outflow), likely in the form of a jet or cone. The GHANG-ES observations 
serendipitously overlap in time with the hard X-ray light curve obtained by Nikolajuk & Walter 
(2013) which they interpret as due to a tidal disruption event (TDE) of a super-Jupiter mass 
object around a 10® Mq black hole. We outline a standard jet model, provide an explanation for 
the observed circular polarization, and quantitatively suggest a link between the peak radio and 
peak X-ray emission via inverse Gompton upscattering of the photons emitted by the relativistic 
electrons. We predict that it should be possible to resolve a young radio jet via VLBI as a result 
of this nearby TDE. 

Subject headings: galaxies: individual (NGC 4845) — galaxies: active — galaxies: jets — galaxies: nuclei 


^Dept. of Physics, Engineering Physics, &; Astronomy, 
Queen’s University, Kingston, Ontario, Canada, K7L 3N6 
irwinSastro.queensu.ca, henriksnQastro.queensu.ca, 
twiegertOastro.queensu.ca . 

^Max-Planck-Institut fiir Radioastronomie, 

Auf dem Hiigel 69, 53121, Bonn, Germany, 
mkrauseOmpifr-bonn.mpg.de. 

^Dept. of Astronomy, University of Massachusetts, 
710 North Pleasant St., Amherst, MA, 01003, USA, 
wqdOastro.umass.edu. 

^US Planck Data Center, The California Institute 
of Technology, MC 220-6, Pasadena, CA, 91125, USA, 
emurphyOipac.caltech.edu. 

^Netherlands Institute for Radio Astronomy (AS- 
TRON), Postbus 2, 7990 AA, Dwingeloo, The Netherlands, 
healdOastron.nl. 

^Physics and Space Sciences Dept., Florida Institute of 
Technology, 150 West University Boulevard, Melbourne, 
FL, 32901, USA, eperlman@fit.edu. 


1. Introduction 


The discovery of a hard X-ray source at the 
center of the galaxy, NGC 4845, by INTEGRAL 
(IGRJ12580+0134) has been interpreted as the 
tidal disruption of a super-Jupiter by a massive 
black hole ( [Nikolajuk fc Walter 20131. As part 
of the Gontinuum Halos in Nearby Galaxies - 
an EVLA[^ Survey (GHANG-ES), we have de¬ 
tected a variable radio source (a compact core) 
in NGG4845 (Table [^, showing unambiguously 
that this galaxy harbours an active galactic nu¬ 
cleus (AGN). The peak of the X-ray light curve 
occurred on January 22, 2011. Our radio obser¬ 
vations were carried out approximately one year 


^The Expanded Very Large Array is now known as the Jan¬ 
sky Very Large Array. 
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later (Table and overlap with the time period 
of their light curve (their Fig. 8 ). 


Table 1: Galaxy parameters 


Parameter 

Value 

Distance (Mpc)*^ 

17 

Inclination (deg)'’’ 

81 

Vq (km s”^)'’ 

1098 

A V (km s”'^)'' 

374 ± 0.4 

HI mass (Al©)® 

(2.1 ± 0.4) X 10® 

Total mass {MqY 

9.9 X 10i° 

Ico (K km s“^)® 

6.07 


° Solanes et al 

|l|2002 

^ Irwin et al. 

2012a 


'^Heliocentric velocity jSpringob et al. |2005| |. Note that 
this (and subsequent HI related parameters) are 
adjustments from [Irwin et al~| l |2012^ given the newer 
HI reference. 


"^From the double-horned profile of Springob et al. (2005i, 
corrected for redshift stretch, instrumental effects and 
smoothing. 


^Springob et al. | | |2005| |, as corrected for pointing offsets, 
source extent, and HI self-absorption. 


•^Using AV and the inclination, above, assuming a spherical 
mass distribution, and an HI radius equivalent to the 
optical blue diameter of 4.8 arcmin (the HI angular extent 
is unknown). 


^CO integrated intensity using Five Colleges Radio 
Astronomy Observatory data ( [Kornugi et al. |2008[ l 


The radio spectrum peaks at Gigahertz fre¬ 
quencies, like the well-known Gigahertz-peaked 
spectrum (GPS) sources seen at high redshift (e.g. 
Tornikoski et al. HMTjlOUeallTQisl ). Character¬ 


istic of such sources is their small radio sizes, gen¬ 
erally thought to be because they are young radio 
sources, or confined by a surrounding medium, or 
both ( Fanti |[2009 Mingaliev et al. ||2013 |. 


NGC 4845 is located in the ‘Virgo Southern Ex¬ 
tension’ (Tully 1982 Karachentsev & Nasonova 


20131 (D = 17 Mpc, fable [^, providing a unique 


opportunity to study the radio/X-ray connection 


for this nearby AGN. A similar relevant case (Za- 


uderer et al. 2011) is the radio/X-ray analy¬ 


sis discussing the ’birth of a relativistic outflow’ 
in a 7 -ray transient at a redshift oi z = 0.354. 
NGC 4845 is 100 times closer and, with future 
monitoring, may provide an unprecedented oppor¬ 
tunity to follow the development of an outflow in 


what could be a present day lower luminosity ver¬ 
sion of a newly forming quasar. 

In this paper, we present our observations and 
develop a simple cone/jet model that is consistent 
with the data, including the variability of the radio 
flux density and spectral index, the behaviour of 
the radio spectrum, the observed circular polariza¬ 
tion, and the absence of significant linear polariza¬ 
tion. We also suggest a possible link between the 
peak X-ray and peak radio emission in the context 
of our model. 


Table provides imaging parameters from our 
five data sets as well as a summary of previous 
radio observations from the literature, clearly il¬ 
lustrating the variability of the source. It is not 
surprising that NGC 4845 has been variously listed 
as having an HII/LINER type nucleus (provided 
by the NASA extragalactic database (NED) based 
on optical spectra of [Ho, Filippenko, fc Sargent 


1995), a pure HII region type nucleus (Filho et al. 
2000|, and a Seyfert 2 nucleus (Veron-Cetty & 


Veron 2010), given that the relative importance 


of the nucleus/disk changes with time. In that re¬ 
spect, NGC 4845 appears to be a closer, fainter 
version of the recently identified ‘changing-look 
quasar’ in which the nuclear activity has changed 
from Type I to a Type 1.9 AGN over a 10 year 
time span ( [LaMassa et al. ||2015 ). 

Sect. presents the observations and data re¬ 
ductions, Sect. |3.1 1 outlines the total intensity re¬ 
sults, in-band spectral indices, and variability of 
the source with fitted spectra, and Sect. ^ gives 
the polarization results. Sect. [^discusses the core. 
We provide a detailed analysis of the core spec¬ 
tra and their variability and also present several 
appendices (Appendices and |H which, not 

only form a description and are consistent with 
the observations, but also summarize more gen¬ 
eral concepts that may be useful for future obser¬ 
vations of NGG4845 and other such sources, es¬ 
pecially in an era in which the sensitivity, fidelity 
and resolution of radio data are rapidly improving. 
The conclusions are presented in Sect. [5] 


2. Observations and Data Reductions 

As part of the GHANG-ES program, five ob¬ 
servations were carried out (B, G, and D configu¬ 
rations at L-band and G and D configurations at 
G-band) in all polarization products over approx- 
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Table 2 

Image parameters 


Observation 

Date^ 


Beam parameters'^ 

Pixel size 

SC Iterations'^ 


rms® 








I 

Q,u 

V 



(GHz) 

//, //, deg. 

// 


{l-iJy beam 


CHANG-ES 

D array L-band 

30-Dec-2011 

1.57470 

38.58, 34.27, -5.22 

5.0 

2A&:P 

40 

27 

28 

C array L-band 

30-Mar-2012 

1.57484 

12.18, 11.10, -41.76 

1.0 

2A&P 

45 

19 

26 

B array L-band 

ll-Jun-2012 

1.57499 

3.51, 3.33, 22.69 

0.5 

2A&P 

18 

15 

15 

D array C-band 

19-Dec-2011 

5.99833 

10.98, 9.06, -1.40 

1.0 

6P+TASP 

15 

15 

12 

C array C-band^ 

23&25-Feb-2012 

5.99854 

3.05, 2.75, -11.71 

0.5 

1P-I-5A&P 

3.9 

3.2 

3.3 

Previous values 

D array Lband (NVSS)^ 

27-Feb-1995 

1.4 

45, 45, 0 

15.0 


450 



Arecibo^ 

2S-Aug-1975 -)■ ll-Sep-1975 
21-Dec-1975 -)■ 03-Jan-1976 

2.380 

~ 162‘ 






B array Lband (FIRST)^ 

09-Oct-1998 

1.4351 

6.4, 5.4, 0.0 

1.8 


150 



C array X-band^ 

21&25Aug-1997 

8.4 

7.08, 5.69, 64.51 



100 



Molonglo^ 

Feb ->■ Aug 1978 

0.408 

172, 172, 0 







Note. —New CHANG-ES data are separated from prevously measured values (where available) by the double-line. 

^Date (or range) that the observations were carried out (UTC, if time specified). 

*^Central frequency of each image. Slight differences for the same band result from different amounts of flagging in different data sets. The 
bandwidths (GHz) were 1.247 1.503 and 1.647 —>■ 1.903 at L and C-bands, respectively. 

^Major Sz minor axis diameters and position angle of the synthesized beam. 

'^Each application of a self-calibration iteration acted on the non-self-calibrated data, but using a model that improved with each iteration. 

®rms map noise measured in a consistent fashion for all maps from random regions that do not appear to contain background sources. For D/L 
and C/L total intensity images only, the rms noise is higher by a factor of 1.5 to 2 close to the source itself because of residual sidelobes and the 
presence of many background sources. 

^These observations were split over two dates as specified. 

^Downloaded data from the NRAO VLA Sky Survey (NVSS) | |Condon et al. |1998[ | website (http://www.cv.nrao.edu/nvss/) 

^Dressel Sz Condon (1978). The observations were carried out during one of the two date intervals specified. 

‘Assuming a zenith angle of 10 degrees. 

■'Downloaded data from the Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey ijBecker et ah | |l995|l website 
(http://sundog. stsci.edu/index.html). 

^Filho et al. | | |2000[ l, using beam parameters corresponding to their most reliable flux density. Calibration errors are quoted as a few percent. 

|Harnett~~| l |l982[ |. 
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imately a six month period (Table [^. We will 
refer to each observation according to ‘configura¬ 
tion/band’ (i.e. B/L, C/L, D/L, C/C, and D/C). 
Data reduction and imaging details are described 


Irwin et al. (2013) and Wiegert et al. 


Here we provide a brief outline as well as details 
specific to NGC 4845. The Common Astronomy 


Software Applications (CASA) package (McMullin 
et al. |[2007 ) was used throughou10 


For each observation, the galaxy was observed 
within a scheduling block (SB) that included other 
galaxies so that the uv coverage could be maxi¬ 
mized for each. For C/C, however, the galaxy was 
observed in two separate SBs on different dates 
which were calibrated separately and then com¬ 
bined for imaging. 

The calibrators were 3C 286 (the ‘primary’ 
calibrator), also used for the bandpass and po¬ 
larization angle calibration, and J1246-0730 (or 
J1224-I-0330 for C/C) (the ‘secondary calibrator) 
which is near the source on the sky. To determine 
the leakage terms between the R and L circularly 
polarized feeds for polarization imaging, we used 
the unpolarized source, J1407-I-2827 = OQ 208, 
except for C-band for which we used the secondary 
calibrators which spanned a parallactic angle > 
60° over the observation^ 


At L-band, the frequency coverage was (in 
GHz) 1.247 ^ 1.503 and 1.647 ^ 1.903 (500 MHz 
total), and at C-band, it was 4.979 —>■ 7.021 (2 
GHz). The frequency gap at L-band was set to 
avoid very strong persistent interference in that 
frequency range. 

Flagging was carried out iteratively, i.e. af¬ 
ter each round of calibration, the data were re¬ 
inspected, more flagging was carried out if needed, 
the calibration tables remade, the data then re¬ 
calibrated, etc. Any flagging that was carried out 
for the total intensity data (RR and LL circular 
polarizations) was also automatically applied to 
the cross-hands (RL and LR). In addition, for the 
polarization calibration (Q and U maps, see be¬ 
low), RL and LR data were inspected, and flag¬ 
ging/calibration iterations were again carried out 


^see also http://casa.nrao.edu 

®The secondary calibrator was used in the event that OQ 208 
carried some polarization at C-band. A test using both 
methods, however, showed negligible differences between 
the results. 


prior to polarization imaging. 


Wide-field imaging (Cornwell, Golap & Bhat- 
nagar [2008 1 using the ms-mfs (multi-scale, multi¬ 


frequency-synthesis) algorithm ( Gornwelll 2008 


Rau & Cornwell 2011) was carried out, with si¬ 


multaneous fitting of a simple power law (S^ oc 
ly'^) across the band. Since NGC 4845 is point-like 
(Sect. 3.1), however, the use of many scales was 
not warranted. After a variety of trials, we found 
that a classic clean gave the best results for highly 
unresolved cases (D/L, C/L, D/C) with the ad¬ 
dition of modest scales for slightly resolved cases 
(B/L: 0, 5, 10, and 20 arcsec; C/C: 0 and 5 arcsec). 

For all maps, Briggs robust = 0 uv weighting 
was used (Briggs 19951, as implemented in CASA. 
A variety of self-calibration steps were then em¬ 
ployed, either phase-only (P) or amplitude and 
phase together (A&P) depending on the specific 
data set and what was required to improve the 
dynamic range of the map. Care was taken to 
ensure that, with successive iterations, the self¬ 
calibration did not go so deep as to start includ¬ 
ing spurious values such as negatives or residual 
sidelobes and the peak specific intensity did not 
decline. The primary beam (PB) correction was 
then applied. Table [^summarizes the imaging pa¬ 
rameters. 


For each observation, imaging of Stokes Q, 
U, and V was also carried out for each observa¬ 
tion using the same self-calibration tables and the 
same input parameters as for total intensity (ex¬ 
cept B/L for which the scales were reduced to 0, 
5, and 20 arcsec). Maps of linear polarization, 
F = were then made, corrected for 

the bias introduced by the fact that P images 


do not obey Gaussian statistics (e.g. Simmons fc 
Stewart [1985 Vaillancourt [2006 1 as well as maps 


of polarization angle. P and V maps were then 
corrected for the primary beam. 


In the case of circular polarization, we also re¬ 
made the maps without fitting a spectrum across 
the band, essentially averaging all spectral chan¬ 
nels together without attempting to fit a spec¬ 
tral index. However, those maps were significantly 
worse than when a spectral fit was carried out. For 
example, at G/L, including a spectral fit increased 
the peak specific intensity by 5% and decreased 
the rms noise by 32%. We therefore retained the 
V maps that included a spectral fit. 
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For each observation, in-hand spectral index, a, 
maps and a error maps were also made for each 


in the CASA clean algorithm. Following this, the 


post-imaging prescription of Wiegert et al. (2015) 


was followed which ensures that the images are 
corrected for the spectral index that is imposed 
by the primary beam (e.g. see also Irwin et al. 


2012b 2013) and smooths over artifacts. A de¬ 


parture from that prescription, however, is that 
we have been more conservative in our cut off, i.e. 
spectral index maps for NGC 4845 have been cut 
off wherever the total intensity falls below lOcr in¬ 
stead of 5 ct (the latter is the norm for CHANG¬ 
ES). 

A weighted mean spectral index, a, and its as¬ 
sociated error, Aa, were calculated over a region 
of interest, taken to be the FWHM of a single beam 
only so, essentially, we are adopting a measure of 
a in a single beam at the galaxy center, a region 
in which variations in the spectral index have been 
minimized. The values are calculated via. 


ite:) 


Aa = 


N, 


1/2 


p/b 


(w) 


( 1 ) 


( 2 ) 


Rau & Gornwell 

2011 

), as implemented ( 

Wiegert et al. |2015 


shown that the true error in the spectral index 
for a typical GHANG-ES galaxy is « 20% higher 


data sets, if we shift the measuring region up or 
down slightly from the center, then a will change 
because of the north-south gradient observed for 
those two data sets (Fig. [^. We have ensured 
that all data sets are measured at the same po¬ 
sition and in the same way and in Sect. |3.1.5| we 
outline a number of checks indicating that possible 
variations in alpha do not affect the conclusions of 
this paper. 

Spectral index maps for the circularly polarized 
emission (not shown) were formed and measured 
in the same way. 

3. Results 

3.1. Total Intensity Emission 

In Fig. we show simple contours (lOcr and 
half-maximum only) of the L-band emission for 
all 3 array configurations on an r-band image of 
NGG 4845 from the Sloan Digital Sky Survey 
(SDSS)j^ It is clear that the radio emission is unre¬ 
solved (cf. the beams at lower right) with the ex¬ 
ception of a small disk around the nucleus showing 
up at the highest resolution (B/L). We will refer 
to this disk as the ‘central disk’. In no case, how¬ 
ever, do we see emission from an extended disk 


(Sect. 3.1.6), i.e. over a size that might approach 


the optical extent of the galaxy. 


where the summations are over the map pixels 
in the region of interest (i.e. the central beam). 
The factor, Np is the number of pixels per beam 
which increases the error since the pixels are not 
independent. These weightings are adopted to en¬ 
sure that pixels with the highest S/N and lowest 
uncertainty are given the highest weight. Since the 
source is mainly unresolved (Sect. [3T| ), the weight¬ 
ing is close to that of a Gaussian. The beam size 
is shown at lower left in Figs. [^and|^ 

The results are given in Table It is impor¬ 
tant to point out that the dynamic range (DR) of 
each total intensity map is very high. The DRs are 
5600/1 (D/L), 5460/1 (G/L) and 11,600/1 (B/L) 
at L-band (but note that for D/L and G/L, the rms 
is higher close to the source, see notes to Table |^. 
At G-band, the DR is higher still, namely 28,300/1 
(D/G) and 88,900/1 (G/G). As a result, the for¬ 
mal errors are quite small. However, tests have 


Fig. i shows G-band contours on a blow-up 
of the same optical emission. This image clearly 
shows asymmetric optical emission in the shape of 
a cone to the north of the nucleus. The galaxy is 
slightly inclined (Table such that its southern 
edge is closest to us; consequently, a correspond¬ 
ing southern optical cone, if it exists, would be 
obscured by the disk. 

Again, the D/G (red contours) radio emission is 
unresolved but the C/G emission (black contours) 
delineates the small central disk quite well. The 
disk subtends 21.7 arcsec (1.8 kpc in diameter) to 
the displayed lOcr contour level at a position angle 
of 83°, in contrast to the larger-scale optical disk 
which shows some distortion and whose position 
angle varies from 74 to 80 degrees, depending on 
the isophotes to which it is measured. A similarly 


■^http: //www.sdss.org/ 
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Table 3 

In-band Spectral Indices 




Date 




Ubservation 

OL 

OLAGN 

OlY 

D 

array L-band 

30-Dec-2011 (Tl) 

+0.95 ± 0.01 

+ 1.11 ± 0.02 

-2.2 ± 0.1 

c 

array L-band 

30-Mar-2012 

+0.97 ± 0.01 

+ 1.11 ± 0.02 

-2.5 ± 0.2 

B 

array L-band 

ll-Jun-2012 

+0.81 ± 0.17 

+0.81 ± 0.17 

-3.4 ± 0.1 

D 

C 

array C-band 
array C-band 

19-Dec-2011 

23&:25-Feb-2012 

-0.453 ± 0.003 
-0.493 ± 0.001 

-0.448 ± 0.006 
-0.493 ± 0.001 



Note. —Weighted mean spectral indices of the total intensity emission, q, the unresolved 
core only, ckagn , and the circular polarization, ov- Measurements were made only for the 
region covered by the half-power beam width at the galaxy center. For quoted error for a 
is the formal random error only; typical errors are likely ~ 20% higher (see Sect.j^. The 
quoted error for OiAGN has been doubled from the pr evious column for resolutions which 
required an adjustment from q to qagn (Sect. |3.1.3|l. For cW, the uncertainties may be 
much higher than the formal error and thus all values of are taken to be roughly the 
same. 


angled radio disk is visible at 8.4 GHz, whose size 
(« 20 arcsec) agrees with our result (Filho et al. 


20001 . 


The central disk shows some extensions to the 
north and south (B/L and C/C, F igs. [T] and 
which could be related to outflow. IFilho et al. I 
(2000) similarly note a small northern extension 
at 8.4 GHz. 


3.1.1. Flux Densities 

Flux densities were measured from the primary 
beam corrected images and are given in Table 
along with previous literature values. With the 
exception of the D/C data set (see below), the 
quoted errors are dominated by variations in flux 
density that result from adjusting the box size 
around the source. It is abundantly clear that 
NGC 4845 is a variable radio source, its L-band 
flux density, S having increased approximately 6- 
fold since the NVSS survey results in 1995. 

Since we will also argue that the flux density has 
varied over the several-month timescale of our own 
observations, it is important to consider possible 
uncertainties that might not have been taken into 
account from the map measurements. These in¬ 
clude a) missing flux due to missing spatial scales, 
b) calibration errors, and c) changes in flux density 
from self-calibration. 

For a), the largest angular scale detected at L- 
band is 16 arcmin and, at C-band, is 4 arcmin. 
The structures shown in Figs. and are well 
within these limits, indicating that all flux has 


been detected in each configuration/frequency- 
band combination. Calibration errors (b), as mea¬ 
sured during the transfer from the primary to 
the secondary calibrator, are less than the errors 
quoted in Table except for D/C for which they 
dominate and are quoted in the table. The un¬ 
certainties of Table |4] do not take into account 
uncertainties in the primary calibrator (3C 286) 
itself (Perley & Butler 12013 of order one per¬ 


cent), but this will not affect comparisons between 
data sets. Self-calibration (c) can result in some 
modification of the flux density. The largest such 
change occurred for the C/L data set, for which a 
4 mJy difference is observed between the non-self- 
calibated data (256 mJy) and the self-calibrated 
data (260 mJy), an increase of 1.5%. Only a de¬ 
crease in flux density would be a concern and any 
decrease observed had a smaller relative error than 
this. Self-calibration, for all data sets, resulted 
in an increase in dynamic range and a lower rms 
noise. 


Since the emission was point-like, the flux den¬ 
sities were also measured by fitting a Gaussian 
to the emission; the result agreed with the tab¬ 
ulated values, except for the high resolution data 
(B/L and C/C) for which some disk-emission is 
also seen and therefore the Gaussian fitted flux 
densities slightly underestimate the total flux den¬ 
sity (92% and 98%, for B/L and C/C, respectively. 
This also allows us to separate the flux density of 
the non-varying disk from the AGN, thereby iso¬ 
lating the flux density of the variable AGN alone. 
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Table 4 
Flux densities 


Observation 

Frequency 

(GHz) 

Time^ 

(days) 

s'’ 

(mJy) 

Sagn^ 

(mJy) 

pd 

(mJy) 

(mJy) 

p/s' 

% 

V/S' 

% 

D array L-band 

1.57 

0 

230 ± 2 

211 ± 3 

0.55 ± 0.05 

6.6 ± 0.2 

0.24 ± 0.02 

2.9 ± 0.1 

C array L-band 

1.57 

91 

260 ± 2 

241 ± 3 

0.35 ± 0.03 

5.7 ± 0.1 

0.13 ± 0.01 

2.2 ± 0.06 

B array L-band 

1.57 

164 

238 ± 1 

219 ± 4 

1.2 ± 0.5 

5.5 ± 0.4 

0.5 ± 0.2 

2.3 ± 0.2 

D array C-band 

6.00 

-11 

432 ± 2 

425 ± 3 

2.3 ± 0.9 

< 0.036 

0.5 ± 0.2 


C array C-band 

6.00 

56 

362 ± 1 

355 ± 2 

0.4 ± 0.2 

< 0.010 

0.11 ± 0.06 


D array L-band (NVSS)® 

B array L-band (FIRST)® 
Arecibo^ 

C array X-band^ 

Molonglo*’ 

1.4 

1.4 

2.38 

8.4 

0.408 


46.0 ± 1 
33.9 ± 0.4 
31 ± 3 
12.5 

193 ± 19 







Note. —Flux densities for the data sets of Tablef^ Frequencies have been rounded. Upper limits assume that the source is unresolved 
and that a detection would require a 3<7 signal, vmere cr is the rms noise of the map (Table [^. 

^Time in days since T1 — 30-Dec-2011. 

*^Flux density of the total intensity emission. The quoted error does not include the 1% absolute flux density error of 3C 286. See 
Sect, f^l.ll for details. 

^Flux density of the AGN after subtracting off the disk e mission . The latter was determined by subtracting off Gaussian fits to the 
peak of the highest resolution data at each frequency (Sect. [3. 1. l[ l. 

Linearly polarized flux density, P = ^, where P has been corrected for positive bias. P has been measured in varying- 
sized regions and, where the emission was point-like (D/L and C/L), via Gaussian fits as well. The quoted error reflects the resulting 
variations. 

®Flux density of the circularly po lariz ed emission. Errors are random only, i.e. as evaluated from the rms noise and residual 
sidelobes in the map itself. See Sect. [3^ 

^Percentage polarization of either the linearly or circularly polarized emission. For V/S, the C-band results have been omitted since 
the V values are non-detections. 

^Flux densities from downloaded data sets (FIRST and NVSS) were measured i n the same way as fo r our data. The NVSS value 
agrees with the result from Condon et al. (1998) and the FIRST value agrees with [Becker et al. [ ) [2012[ l. 

^See Table for references. 
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Fig. 1.— L-band contours on a Sloan Digital Sky Survey (dr7) r-band image of NGC4845. Contours are set 
at lOtT and 50% of the peak for each data set. D/L: grey (0.40, 112.4 mJy beam“^), C/L: black (0.45, 122.8 
mJy beam“^), and B/L: blue (0.18, 104.7 mJy beam“^). The beam sizes are shown at lower right. 
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Fig. 2.— C-band contours on a blow-up of the optical emission as shown in Fig. Contours are red for 
D/C at 0.15 (lOcr) and 212.4 mJy beam“^ (half-peak value). Contours are black for C/C at 0.039, 0.10, 
0.30, 1.00, 6.00, and 173.1 mJy beam“^ (half-peak value). The beam sizes are shown at lower right. 
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Fig. 3.— L-band spectral index maps formed as described in Sect. Colours represent the spectral index 
with the colour legend shown at right. Grey solid contours represent values from the corresponding error 
map, with smallest values in the interior to larger values towards the outside. Note that the size and colour 
scales are different from map to map. The ellipse at lower left shows the beam. Left: D/L spectral index 
map with error contours at 0.01, 0.025, and 0.05. Center: C/L spectral index with error contours at 0.01, 
0.025, and 0.05. Right: B/L with error contours at 0.15 and 0.2 (south-west contour is 0.15 and broken 
north contour is 0.2. 



J2000 Right Ascension 


Fig. 4.— C-band spectral index maps formed as described in Sect. Colours represent the spectral index 
with the colour legend shown at right. Grey solid contours represent values from the corresponding error 
map, with smallest values in the interior to larger values towards the outside. Note that the size and colour 
scales are different from map to map. The ellipse at lower left shows the beam. Left: D/C spectral index 
map with error contours at 0.01, 0.05, and 0.1. Right: C/C with error contours at 0.01, 0.1, and 0.2. 


Sagn (Table I^. From the global flux densities of 
the disk alone at the two frequencies (19 ± 4 mJy 
at L-band and 7 ± 1 mJy at C-band) we find that 
the spectral index for the disk is = —0.74. 

Such a value is typical of galaxy disks, in general. 

3.1.2. Core Size and Astrometry 

The Gaussian fits discussed in the previous sec¬ 
tion were also used to determine the size and as¬ 
trometry of the central core (which we take to rep¬ 
resent the AGN), using the highest resolution data 


(C/C and B/L, Table [^. The average central po¬ 
sition of the source is at RA = 12*'58™01.1968® ± 
0.0007" (±0.01"), DEC = 01°34'32.42" ± 0.01", 
where we adopt the higher of the uncertainties 
from the two data sets. This is the most precise 
position yet determined for the core of NGC 4845 
and should also be the most accurate, given that 
radio data do not suffer from extinction as do op¬ 
tical data. Nevertheless, our position agrees well 
with that of the SDSS (uncertainty of 0.2") as 


®These are formal errors from the CASA imfit task. 
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Table 5: Source astrometry 


Observation 

(GHz) 

RA, DEC (J2000)“ 
h:m:s, o:/:// 

Source size^ 

//, deg. 


B array L-band 

1.57499 

12:58:01.1955 ± 0.0007, 4-01:34:32.42 ± 0.01 

0.67 ± 0.03, 0.15 ± 0.09, 87 ± 

170 

C array C-band 

5.99854 

12:58:01.1981 ± 0.0001, 4-01:34:32.424 ± 0.001 

0.40 ± 0.03, 0.36 ± 0.03, 97 ± 

175 


“Reference frequency for each image, as in Table |2] 

^'Positions from Gaussian fits. The quoted uncertainty is the larger of the formal errors or variations resulting from the adoption 
of three different boxes within which the fit was carried out. 

“Deconvolved from the beam, from Gaussian fits. The quoted uncertainty is as described in note b. The fit includes an offset to 
account for the fact that the AGN is within a disk. 


quoted in NED. 

The size of the central core, after deconvolving 
from the synthesized beam, is also provided. A 
comparison with the beam sizes given in Table 
shows that the central core is unresolved, as mea¬ 
sured in the highest resolution data sets. (Repeat¬ 
ing this exercise for the lower resolution data gives 
the same conclusion.) The deconvolved source size 
(0.2 —0.7 arcsec) is also of order or smaller than 
the pixel size (0.5 arcsec = 41 pc) which sets a 
limit on the possible precision of size measurable 
by a Gaussian fit. Therefore, the unresolved core 
could be smaller still. 


3.1.3. In-band Spectral Indices 


Since the source is variable (Sect. 3.1.41, a ‘clas¬ 


sical’ computation of spectral index between L- 
band and C-band would not yield reliable results 
since the two different bands were measured at dif¬ 
ferent times. Fortunately, the wide bands used in 
the CHANG-ES observations make it possible to 
compute in-feandspectral indices (Sect.[^. Figs.[^ 
and 1^ show the resulting maps for L-band and C- 
band, respectively. 

Fig. shows successively higher resolution (left 
to right) L-band spectral index maps. A largely 
north-south gradient is visible in the lower reso¬ 
lution D/L and G/L images, though not in the 
high resolution B/L data. These gradients are in¬ 
triguing, since they appear to align in the direc¬ 


tion of the optical cone (Sect. 3.1) and since they 
are visible in two independent data sets. Since 
the gradients become most pronounced after self¬ 
calibration, a variety of tests were run with dif¬ 
ferent self-calibration iterations, inputs and box 
sizes; nevertheless, the gradient remains. To be 
cautious, however, we have computed the weighted 


average spectral indices, a (see details. Sect. [^, 
only for the central FWHM of the emission shown 
in Fig.|^ That is, a applies to a size equivalent to 
the beam size shown at lower left rather than to 
the entire region shown in the figures. This also 
ensures that the spectral indices are measured only 
in regions where the error (as shown by contours) 
is lowest and emission is highest. The results are 
given in Table Note that L-band spectral indices 
in all configurations are positive, and of order -hi. 

Fig.|4] shows the G-band spectral indices. The 
weighted mean spectral indices, again for the cen¬ 
tral FWHM of the emission, are given in Table 
C-band spectral indices are all negative, of order 
-0.5. 


It is clear that the spectral index of this com¬ 
pact core peaks between L-band {v — 1.57 GHz) 
and G-band {v = 6.00 GHz). The source there¬ 
fore may be a nearby lower-luminosity analogue 
of the GPS sources seen at higher redshifts (e.g. 
O’dea, Baum, fc Stanghellini ||1991 ). 

The importance of measuring in-band spectral 
indices cannot be overstated for these observa¬ 
tions. Since the source is variable, a classical mea¬ 
surement of spectral index between the centers of 
the two bands would have given grossly erroneous 
results, even if the observations had been simul¬ 
taneous, since a straight line (in log space) would 
simply have been drawn between the flux densities 
of the mid-points of the two bands (see also Fig.|^. 
Our in-band spectral indices correspond to the 
slopes of the spectrum at the midpoints of both 
bands, a fact that we will exploit in Sect. |3T75 


Finally, the spectral index of the weak disk 
emission, a^isfc = —0.74 (Sect. 3.1.1) is not the 
same as that of the central core in either band. 
Since the HPBW region in which a was measured 
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for the low resolution data (D/L, D/C, and C/L) 
encompasses both the core and the disk, a small 
correction is required to obtain the spectral index 
of the core itself. For each band, this was done 
by determining the flux density of the disk at the 
upper and lower ends of the band using adisk, de¬ 
termining the total flux density for the upper and 
lower ends of the band using o, subtracting off the 
disk flux density from the total at the two band 
ends, and recalculating the in-band spectral in¬ 
dices, yielding uagn- These results are also given 
in Table |3] 


3 . 1 . 4 . Variability 

Table reveals a significant change in flux den¬ 
sity with time for NGC 4845. 

For the historical data at L-band, the flux den¬ 
sity appears to decrease from the 1995 NVSS (D- 
array) value of 46 mJy to the 1998 FIRST (B- 
array) value of 34 mJy. However, the FIRST B- 
array data reveals only a point source without the 
disk since the rms noise is of order the brightness 
of the disk as measured in the CHANG-ES B/L 
data set. The difference between those two flux 
densities (12 mJy) may therefore be attributable 
to some missing low brightness extended emission 
of the central disk which we have measured to be 
19 ± 4 mJy (Sect. [3TTT ). Two additional earlier 
archival VLA data sets are also available, though 
unpublished, which we have now reduced. These 
are L-band D-array data (Mar. 23, 1988) and 
L-band B-array data (Nov. 15, 1987) for which 
the flux densities are 48 mJy and 28 mJy, respec¬ 
tively. Given uncertainties (especially at the ear¬ 
lier epoch), we consider these values to agree with 
the NVSS and FIRST results, respectively. There¬ 
fore, of the limited data available, there was no 
detected L-band variability prior to 1998. 

On the other hand, our GHANG-ES observa¬ 
tions have now shown an increase at L-band of 
approximately a factor of 6 compared to the NVSS 
results. Unfortunately, there is no prior measure¬ 
ment of G-band flux density in the literature. The 
(1986 to 1987) Green Bank 4.85 GHz survey (GB 6 , 


Gregory et al. 19961, which covered this region 
of sky to a limit of « 42 mJy (their Fig. 2), 
did not detect the source. Allowing for the 5 ct re¬ 
quirement for the source to be included in their list 
suggests that the G-band flux density was < 210 
mJy at the earlier epoch, or a factor of at least 2 


lower than our peak G-band measurement. We 
suspect that the flux density change was much 
higher, though, given the L-band measurements 
and the fact that the source is highly absorbed at 
L-band but not G-band (Sect. [T^ . 

A plot of AGN flux density and in-band spec¬ 
tral indices (both corrected for the disk) as a func¬ 
tion of time for the GHANG-ES data is shown in 
Fig[^ Earlier historical data have been omitted, 
given the large lapse in time. This plot shows 
that the AGN flux density has varied even over 
the few month time period of our observations. 
At L-band, the flux density increases, and then 
decreases. At G-band, the flux density shows 
a strong decline. The spectral indices (dashed 
curves) also reveal variations in both bands. 

Erom these data, we adopt 3 time stamps, 
namely Time 1 (Tl): the first time at which L- 
band data are obtained and for which G-band data 
can be interpolated from its curve. Time 2 (T2 = 
Tl -I- 56 days): the last time at which G-band 
data have been obtained and for which L-band 
data can be interpolated, and Time 3 (T3 = Tl -f 
196 days): an extrapolation of both L-band and 
G-band data to the time at which the flux densities 
at both bands are equal. In all cases, simple lin¬ 
ear interpolations or extrapolations are performed, 
and the same is done for the spectral indices. The 
resulting data (Table provide flux densities and 
spectral indices for matching times. 

3.1.5. The AGN Spectrum and its Variability 

Erom the data of Fig. it is possible, not only 
to estimate the spectrum of the AGN, but also to 
see how it varies with time. For each time stamp, 
we now have four data points that can be used 
to constrain the spectrum: the L-band flux den¬ 
sity, the L-band spectral index, the G-band flux 
density, and the G-band spectral index, where the 
spectral indices represent tangent slopes of the 
spectrum at the two frequencies. We then fit a 
polynomial to the spectrum of the form 

Si, = ao + axv + a2i'‘^ + 03 (3) 

where S^, is the flux density at frequency v and 
are constants. The derivative of this equation is 
equated to the derivative of the flux density equa¬ 
tion, Su oc v°‘ I both of which represent tangent 
slopes to the spectrum at frequency, v. With two 
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Fig. 5.— AGN flux densities, Sagn (solid lines with circular markers), and in-band spectral indices, a 
(dashed lines with diamond markers), are shown for the L-band (blue) and C-band (green) data, from 
Tables 1^ and 1^ as a function of time, where T1 = 30-Dec-2011. Error bars are shown but are mostly smaller 
than the marker size. Dotted lines show linear extrapolations of the curves to the date, 13-Jul-2012 (T1 + 
196 days), at which the L-band and C-band flux densities are equal (Table [^. The dashed curve shows an 
extrapolation of the C-band flux density of the form. Sc oc (Eqn. 22). 


flux densities and two slopes at each time stamp, 
Eqn. 1^ was solved for the spectrum at each time 
stamp. 

The results are given in Table and the spec¬ 
tra are plotted in Eig. It is now clear that the 
turnover of the spectral index is shifting to lower 
frequencies with time, and the peak of the spectrum 
also declines with time. The peak frequency shifts 
from 4.9 GHz to 4.0 GHz to 3.2 GHz. 

We have carried out some tests to verify the va¬ 
lidity of this result. Eor example, the same trend 
can be seen in the original data prior to subtract¬ 
ing the disk emission. If we repeat this exercise 
prior to subtracting the disk contribution to both 
flux density and spectral index, then Vm is altered 


by less than 1% for all time stamps, Smax for time 
steps, T1 and T2, increase by less than 2% and 
the peak flux density for time step, T3, increases 
by 15%. Given the error bars in the L-band spec¬ 
tral index, we also repeated the analysis assuming 
that ol has not varied at all (i.e. letting aZ ~ IT 
for each time stamp), then there is a 1% decrease 
in the peak frequency, Vm, and a 5% decrease in 
the peak flux density, Smax, for the T3 curve only. 

Eitting a curve, where possible, instead of simple 
straight lines (e.g. using a second order polyno¬ 
mial for the L-band data of Eig. [^, also does not 
alter the above conclusions. 

The final issue is related to extrapolation/interpolation. 
Note that point T1 is measured at L-band and in- 
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Table 6 

Flux Densities and Spectral Indices for Matching Times 


Date 


L-band 



C-band 



Sagn 

Frequency 

OiAGN 

Sagn 

Frequency 

OLAGN 


mJy 

GHz 


mJy 

GHz 


T1 (30-Dec-2011) 

211 

1.575 

+1.11 

414 

5.998 

-0.455 

T2 (24-Feb-2012 = T1 -|- 56 days) 

229 

1.575 

+1.11 

355 

5.999 

-0.493 

T3 (13-Ju 1-2012 = T1 -|- 196 days) 

209 

1.575 

+0.68 

209 

5.999 

-0.587 


Note.— L-band and C-band fl ux den sities and spectral indices for matching times, interpolated or extrap¬ 
olated from Fig. (see also Sect. [^1.4[l. 


terpolated at C-band; point T2 is measured at 
C-band and interpolated at L-band. Only point 
T3 is an extrapolation in time. If we instead 
make a more modest extrapolation (i.e. T1 -|- 100 
days, or 44 days after the last measured C-band 
point but still an interpolation at L-band), then 
the trend is still consistent with what is shown in 
Fig.[§ again with a declining peak that is shifting 
to lower frequencies with time. This intermediate 
curve (not shown) falls between the T2 and T3 
curves of Fig. 

At time T3, the L-band extrapolation is small, 
both when one considers the elapsed time (32 days 
after the last L-band measurement) as well as the 
modest implied variations in flux density and spec¬ 
tral index. It is the C-band flux density that has 
undergone the largest change as a result of our ex¬ 
trapolation (Fig. [^. We have made the simplest 
possible assumption as to the form of its time de¬ 
cay, i.e. linear, arguing that to impose anything 
more complex (and arbitrary) is not justified by 
the data. 


Howev er, in anticipation of our model results 
(see Sect. 4.3) we can also try & (x de¬ 

pendence (Eqn. on the C-band flux density. 


The results are shown as the black dashed curves 
in Figs. and Note that we restrict this curve 
only to the region from T2 to T3 in which a true 
extrapolation at C-band is needed. In addition, 
the curve’s time dependence applies to the highly 
optically thin regime, whereas there is actually a 
changing optical depth with time, as explained in 
Sect. |4.3| (hence the ‘disjoint’ slopes at T2). This 
new extrapolation increases the C-band T3 flux 
density by only 3% from 209 mJy (Table[^ to 215 
mjy and its affect on the behaviour of the curves 
in Fig. is minor. 


Any other more dramatic extrapolation is both 
unknown and speculative and, as stated above, 
even modest extrapolations still follow the trends 
shown in Fig. In summary, we conclude that 
the behaviour of declining flux densities and fre¬ 
quencies with time is a robust conclusion within 
the constraints of the available data. 

Fitting the spectrum with a polynomial is the 
best and most accurate way of describing the 
source’s variable behaviour with time, as observed. 
It does not assume any physical model, however, 
an issue that we will address in Sect. U) 


3.1.6. The Absence of Large-Scale Disk Emission 

No extended disk emission beyond the 1.8 kpc 
diameter central disk seen in Figs and |(B/L 
and C/C data, respectively) has been detected 
from NGC 4845. This does not appear to be a 
sensitivity or dynamic range problem when com¬ 
pared to other galaxies in the CHANG-ES sam¬ 


ple (Wiegert et al. 2015). In addition, measur¬ 


able emission is seen in NGC4845 throughout the 
large-scale disk in all WISE (Wide-held Infrared 
Survey Explorer; IR) bands as well as in GALEX 
(Galaxy Evolution Explorer; UV) images. Previ¬ 
ous radio continuum observations (Table sim¬ 
ilarly reveal radio emission that is spatially con- 
hned to the central region of the galaxy. A thor¬ 
ough discussion of this issue, for example, a study 
of SE rates and timescales in comparison to the 
radio emission, must await a future paper. We 
do note, though, that truncated gaseous disks and 
galaxy stripping are well-known properties of some 


Virgo Gluster galaxies (e.g. Cayatte et al. 1990 


Kenney fc Koopman |T9^ Kenney, van Gorkom, 

fc Vollmer ||2004 Vollmer et al. p013 Murphy et 
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Fig. 6.— Spectra from polynomial fits (Table for three time stamps: T1 (30-Dec-2011) in red, T2 = T1 
+ 56 days in green, and T3 = T1 + 196 days, in blue. Data (black dots) are from Table and the peak 
of each curve is marked with a small bar. Note that each curve is an explicit solution to a polynomial fit 
with 4 constraints: 2 flux density points and 2 slopes. T3 data have been extrapolated from the measured 
data. The black dashed curve corresponds to the non-linear extrapolation shown by the black dashed curve 
in Fig. Measurement error bars are typically smaller than the points. 


ar~|[2^ . 

3.2. Polarized Emission 

Flux density measurements for linear polariza¬ 
tion, P, circular polarization, V, and their frac¬ 
tions of the total unpolarized flux density, S, are 
given in TableIn all arrays and frequencies, the 
percentage linear polarization is less than or of or¬ 
der 0.5%. Since the polarization calibration can¬ 
not be guaranteed to be better than 0.5% (S. My¬ 
ers, private communication)]^ we cannot claim to 
have detected any linear polarization in any con- 


See also the NRAO polarimetry on-line documents. 


figuration/band combination for NGC 4845. 

On the other hand, there is a clear measure¬ 
ment of circular polarization at L-band, though 
not at C-band. The circular polarization is quite 
strong, with a flux density of 5.5 to 6.6 mJy, (cf. 
the rms noise of 15 to 28 /rJy beam“^; Table 
or a polarization percentage of 2 to 3%. Signifi¬ 
cantly, each independent observation has resulted 
in approximately the same percentage circular po¬ 
larization. 

It is important to consider whether the circu¬ 
larly polarized emission could have been artifi¬ 
cially introduced. For example, any errors will 
affect a V signal more strongly than total inten- 
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Table 7 

Parameters of Polynomial Spectral Fits 


Date 

ao 

a-i 

02 

03 

Urn^ 


Time 1 (T1 = 

30-Dec-2011) 

-100.27 

249.79 

-35.17 

1.305 

4.87 

432.7 

Time 2 (T2 = 

24-Feb-2012) 

-156.93 

338.22 

-65.22 

3.845 

4.03 

398.5 

Time 3 (T3 — 

13-Jul-2012) 

-36.78 

230.71 

-53.02 

3.566 

3.22 

275.4 


Note. — Parameters of the polynomial fit to the spectrum for 3 different time stamps. 
The fit is to Su — ao ai i' <12 + aa where Si, is the flux density (mjy) at 

frequency, ly (GHz), and simultaneously matching the slope of this f unction to the slope 
of the power law spectrum, Su oc u°‘ in each band (see also Sect. |3.1.5f. Input data 
are from Table 

^Frequency of the spectral peak (GHz). 

^Flux density at (iriJy). 


sity since V is measuring a difference between cor¬ 
related signal^ It is well known that the VLA 
suffers from ‘beam squint’ which can result in spu¬ 
rious signals off-axis (Brisken 20031. The center 
of NGC 4845, however, was placed at the field 
center where beam squint is zero. Self-calibration 
can also result in spurious V signals if there are 
other off-axis circularly polarized sources in the 
field (see Homan &: Wardle~||1999 for details of a 
related example). 

Our best counterexample, however, is the B/L 
data for which negligible signal was observed other 
than in the centre of the map. This data set 
was re-imaged with no self-calibration table ap¬ 
plied and (although beset with higher residual 
sidelobes due to the lack of self-calibration) the re¬ 
sulting image still shows a significant V signal. An 
alternative self-calibration method was also em¬ 
ployee!^ again, the result was poorer with higher 
sidelobes, but there was no change in the V flux 
density. Imaging of the primary calibrator, 3C286, 
showed V/I < 0.17% and the polarization calibra¬ 
tor, OQ208, had V/I < 0.08%. 

In summary, the circular polarization appears 
to be of order a few percent in V/S and also 
V/Sagn (Table I^. Moreover, we find the V spec¬ 
tral index, ay, again possible to determine be¬ 
cause of our broad bands (Sect. [2), to be quite 
steep, of order —2 —)• —3.5 (Table Given the 


= {RR — LL)I2, where RR refers to the correlation of 
‘right-right’ -hand polarization and LL refers to ‘left-left’. 
®A ‘T’ type self-calibration in which a single solution for 
both polarizations is obtained. 


above discussion, we do not claim that the small 
variations between arrays are significant. The 
rather striking result of measurable circular po¬ 
larization and its variation with frequency will be 
discussed further in Sect. 14.41 


4. Discussion 


In this section, we seek to explain the following 
observations: a) the turnover in the radio spec¬ 
trum and its evolution with time (both ampli¬ 
tude and frequency shifts), conveniently summa¬ 
rized for this ‘Gigahertz-Peaked Spectrum’ (GPS) 
source in Fig. and b) the detection of circular 
polarization while a believable linearly polarized 
signal is absent. We also attempt to connect the 
X-ray observations of Nikolajuk & Walter (2013) 
with our radio data. 


We describe more generally a jet/cone model 
with basic geometry that lies behind our general 
thinking in Appendix]^ For this section, it is suf¬ 
ficient to imagine a source of relativistic particles 
as illustrated in the cone geometry of Fig. The 
optical cone observed in Fig. [^suggests such a ge¬ 
ometry on larger scales, possibly representing an 
earlier outflow. As we do not know the small-scale 
opening angle, however, we will use ‘jet’ or ’cone’ 
interchangeably when referring to our model. 

We are aided by the fact that the X-ray 
outburst was observed and provides us with a 
timescale. We will refer to the time of the peak 
of the X-ray light curve (22-Jan-2011) as <o (the 
disruption of a mass that apparently triggered 
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Fig. 7.— One side of an outflow, of order 0.1 pc in size, which forms a wide-angle cone iO large) or narrow 
angle jet (small 6 ) of relativistic electrons. The origin of the outflow is at the apex and the helical curve 
represents the azimuthal magnetic field. A suggested line of sight (los), which passes through a sheath of 
material in or around the cone, is shown by the heavy back arrow which makes an angle, t, with the radial 
(along the cone axis) direction. A depolarizing screen (not shown) is outside of the cone along the los. 
Inset: Blow-up of the apex region, of order 1000 AU in size, showing the shocked region in which X-rays are 
generated. 


the X-ray flare occurred 60 to 100 days earlier, 
Nikolajuk &: Walter ||2013l ). Our radio data span 
6 months (Table and our hrst timestamp from 
our spectral fits (Table[^ occurs at T1 = 342 days 
« 1 year after the X-ray peak. 

Several reasonable simplifications can be made. 
First, for a very fast jet, transformations would 
be required between the source and observer’s rest 
frame for parameters such as jet angle, brightness, 
magnetic field, etc. For completeness, these trans¬ 
formations are provided in Appendix]^ however, 
we Hnd that unless the bulk jet speed is ultra- 
relativistic, such transformations are not required 
and will not be applied here. Secondly, our calcu¬ 
lations are order of magnitude estimates to test for 
feasibility and consistency with the observations. 


4.1. Source Size 


From observations, the core of NGC 4845 is un¬ 
resolved, likely less than 0.5 arcsec (Sect. 3.1.2) 
so we do not have sufficient resolution to mea¬ 
sure its size, Og, directly. However, the variabil¬ 


ity itself suggests a limiting size. If the source 
was created with the X-ray flare observed by IN¬ 


TEGRAL in January, 2011 (Nikolajuk & Walter 


2013), then the projected physical size should be 
< cto sin (i) « 10^® sin (i) cm, where t is the incli¬ 
nation of the major source axis to the line of sight 
(Fig.[7l) and to is the elapsed time. For to ~ 1 year 
(cf. dates of Table and a distance of 17 Mpc 
(Sect.[^, this suggests that Og < 4sin(t) mas. If 
i = 0°, then, we are looking down the major axis of 
the cone along the outflow and Og would relate to 
the cross-sectional diameter of the cone; if t = 90°, 
then the cone is aligned with the sky plane and Og 
would relate to a typical cone size in projection. 
In the following, we adopt 6 g ^ 1 mas (~ 0.1 pc 
or « 10^^ cm) as a hducial value. 

Note that the brightness temperature at 20 cm 
is « 3.4 X 10^^ K if the source size is of order 1 
mas. This is still below the limit for catastrophic 
inverse Compton losses. GPS sources at higher 
redshift suggest a range of sizes with the lowest « 
1 pc (Fanti 2009), hence the core of NGC 4845 


extends the array of GPS source sizes an order of 
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magnitude smaller. 


4.2. The Low Frequency Spectral Turnover 

A consistently difficult problem in the past 
has been to understand whether low frequency 
turnovers are caused by thermal absorption or 
synchrotron self-absorption (SSA). In this section 
we explore both possibilities and conclude that 
the dominant cause of the turn-over is likely syn¬ 
chrotron self-absorption. 


We first, however, exclude Razin suppression 
of the low frequency emission as the cause of the 
turnover because of the predicted exponential cut¬ 


off below the Razin frequency (Melrose 1980) 


that is not observed. At least for a homogeneous 
source, the effect also does not give a dependence 
of circular polarization on frequency which in fact 
we detect (Sect. 3.2). 


4-2.1. Synchrotron Self-Absorption 

For synchrotron self-absorption (SSA) as an ex¬ 


planation for a GPS source (see e.g. Artyukh et al. 


2008) one expects the brightness temperature of 


the source to approach the electron kinetic tem¬ 
perature (e.g. Longair ||1994 |. The source bright¬ 
ness temperature in terms of the wavelength. A, 
and flux density, S,^, 


IS 


T 

2k n., 


7 X lO^Aj 
0,(mas)2 " 


K, 


(4) 


where 17^ = (7r/4)0g, 6 s is the source size, and k is 
Boltzmann’s constant. We have taken A in units of 
6 cm, ds in mas and the flux density in mJy. The 
kinetic temperature, Tg, for relativistic electrons 
is given by kT^ = yrrieC^/S. Equating these two 
temperatures gives the limiting self-absorbed flux 
density at a peak frequency for electrons of energy, 
7 , as 


S. 


r /2 

V 3 3 


950 


v^^^6s{masY 

i?±(-2)i/2 


mJy. 


(5) 

where Vg = eRj_/(27r)TOeC is the gyrofrequency, e 
and TOe being electronic properties, and is the 
frequency in units of 5 GHz. From this last expres¬ 
sion we can infer the universal (thermodynamic) 
self-absorbed spectral index, a = 2.5 and declining 
flux densities with decreasing frequency. 


Adopting i3j_(—2) = 1 (see Appendix [e|) and 
9s{mas) = 1 (Sect. 4.1), the observed flux density 
at 6 cm is lower than the flux density given by 
Eqn. indicating that self-absorption is not com¬ 
plete at this frequency. However, at 20 cm, the 
observed flux density is much higher than that 
of Eqn. (by an order of magnitude). In other 
words, the observed flux density at 6 cm implies 
a peak (turnover) frequency below 6 cm whereas 
the observed flux density at 20 cm implies a peak 
frequency well above 20 cm. Thus we expect our 
results summarized in Fig. I^to describe a spectral 
region that is transiting towards self-absorption at 
the low frequencies. 

One of the principal results of these observa¬ 
tions is the measurement of the in-band spectral 
indices, which have been used to define the fits of 
Fig.i If we are to explain these curves as a transi¬ 
tion to synchrotron self-absorption, then we must 
consider the roll-over in spectral index from an op¬ 
tically thin spectral index of a = (l—p)/2 « —0.5 
=> p = 2 (Table (using N{E) = NoE~p for 
the energy distribution function of the electrons) 
towards the index a = 2.5 which is expected when 
the source is fully self-absorbed, as shown above. 
This roll-over is described in ( Pacholczyk [1970 p. 
97). We calculate 


din J(z) 
dlnz ’ 


( 6 ) 


where 

J(z) = z 5/2 — exp (— z-(2+P/2))) , (7) 

and 0 = V jvx- The frequency at which the optical 
depth equals one is vi, which for p = 2 is given 
by vi = 0.707i^cj I'c the critical frequency. The 
result is shown in Fig. for several values of p. 
The reference case for our purposes is the lower 
curve at large z for which p = 2. We see that a 
varies from the optically thin value of —0.5 to a 
moderately optically thick value of a = -1-1.0 (as 
observed at L-band) over roughly a factor three 
in frequency. However because of the steepness of 
the rising curve at z = 1 it is possible to have 
a in the optically thick region ranging from just 
below 0.5 to 1.0 over essentially the same range 
in frequency. This is a reasonable explanation of 
the curves in Fig. [^between 6 and 20 cm (G-band 
and L-band, respectively). Given the index p, we 


17 

















emphasize that the only free parameter is vi in 
describing the spectral shape. To detect the fully 
self-absorbed spectral index, one has to descend 
another factor of two in frequency. A fit by the 
synchrotron model to all the data is included in 
Fig.lH 



Fig. 8 .— The figure shows the rollover in the spec¬ 
tral index during the onset of synchrotron self¬ 
absorption at low frequencies. The frequency is 
measured in units of vi, the frequency at which 
the optical depth is 1 (i.e. z = vjvx). This is 
0.707 times the peak frequency for p = 2. The 
lower curve at large frequency is our reference case 
for which p = 2 {a = —0.5). The higher curves at 
large frequency taken in rising order are for p = 1.5 
and p = 1 respectively. The transition from thin to 
moderately thick is over roughly a factor of three 
in frequency. Another factor of two in frequency 
is required to attain the completely self-absorped 
index. 

We can also estimate the relativistic electron 
density required to produce the turnover. The ab¬ 
sorption coefficient for the synchrotron particles is 
approximately 

Kerv = 1-7 X A^q ^ cm~^, (8) 

i/f 

where the subscript, er, refers to relativistic elec¬ 
trons and i^i is in units of GHz. When p = 2, the 
coefficient No « UerEo- Here tier is the relativis¬ 
tic electron density and Eo = ^orrieC^ is the lower 
energy cutoff of the distribution. If L measures 


the line of sight distance through the source then 
we conclude, by setting the optical depth equal to 
one that 


L = 2 X 10 
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1 


7o(10)S^(-2)^ 


cm 


-2 


(9) 


Here 7o(10) is the lower cutoff in units of ten 
(which we set to 1 , following other authors) and 
we have taken vi = 1.5. Nominally this suggests, 
given our proposed source size ( 10 ^^) that Uer ~ 
200 cm~^. 


If the number of protons is the same, then a 
rough calculation of mass for a cylinder (approx¬ 
imating a cone) of length and diameter, 10 ^^ cm, 
would result in ^ 3 x 10^® g (more if thermal 
gas is entrained). This is 5 to 10% of the mass of 


3 X 10^° < M < 6 X 10^° g estimated by Nikola- 


juk & Walter (20131 to have produced the X-ray 
flare via accretion of a ‘super-Jupiter’ object. In 
the latter case, the authors suggest that 10 % of 
a ‘super-Jupiter’ of mass, 14 - 30 Jupiter masses, 
has produced the X-ray flare. This suggests some 
consistency between the relativistic electron den¬ 
sity required to explain the radio spectral turnover 
and the expected quantity of particles supposed to 
be involved in the X-ray event (Fig. [^. 

For a more thorough investigation of SSA as a 
possible cause of the turn-over, we have also fit the 
T2 timestamp with a pure SSA spectrum, allow¬ 
ing vi, p, and the amplitude to freely vary. For 
the latter, we introduce a scaling factor, C, which 
converts the amplitudes to a flux density, S'(z), for 
comparison to the data, i.e. S'(z) = C J{z). The 
best fit is chosen by minimizing during a least 
squares fit. 


The resulting best fit parameters are given in 
Table and the spectrum, along with its deriva¬ 
tive, are plotted as the red curves in Fig. This 
curve is identical to the generic red curve in Fig.[^ 
when i/i = 1.8 GHz. 


4-2.2. Thermal absorption 

In this section we consider the synchrotron 
emission to be optically thin and the turnover to 
be due to thermal absorption only. 

In the simplest approach, since we know the 
optically thin spectral index, which we take to be 
a « —0.5 (Table [^, the spectral index in the 
optically thick (r = 1 ) regime that would result 
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Table 8 

Parameters of Synchrotron and Thermal Spectral Fits at T2 



P 

yi 

k 

c 

Sl 

Sc 

OiL 

Oic 


Synchrotron^ 

2.0 

1.8 


645 

358 

349 

1.20 

-0.460 

118 

Thermal^ 

1.8 


1.8 

740 

299 

344 

1.05 

-0.300 

249 

Synchrotron -|- ThermaF 

2.0 

1.5 

0.3 

675 

336 

329 

1.16 

-0.439 

117 

Data‘s 





229 

355 

1.11 

-0.493 



Note. — Best fit parameters for the various fits as described in Sects. |4.2| and Fig. 
Subscripts, L and C refer to L-band and C-band, respectively, is a measure of the good¬ 
ness of fit, here defined as x^ — {['S’n(/^t) — Sj:,(data)]^ -1- [Scifit] — Sc{data)]^}/(AS)'^ + 
{(/R) — Oil, {data)]‘^ [oicifit) — o-c (data)]^}/(Aa)^, where we have taken AS = 5 mjy 
and Aa = 0.01 as r epresentative values. Note that C is defined differently, depending on 
the fit (see Sect. [T^ . 

^ Low-frequency turnover due to synchrotron self-absorption only. 

'^Turnover due to thermal absorption only. 

^Turnover due to synchrotron self-absorption plus a foreground thermal screen. 

'^Data are from Table 




Fig. 9.— Best fit synchrotron-self-absorption-only spectrum (red), thermal-absorption-only (green) and a 
synchrotron self-absorption plus a thermally absorbing screen (blue), all as applicable to time, T2. Best 
fit parameters are provided in Table Black dots represent the data. The left plot shows the spectra, 
displaying fits to the flux densities, whereas the right plot shows the slopes, displaying fits to the spectral 
indices. Measurement error bars are typically smaller than the points. 


from thermal absorption can be computed. In Ap¬ 
pendix]^ we show that purely thermal absorption 
can not provide a satisfactory fit to both the low 
frequency and high frequency spectral indices. If 


one improves the fit at one end by adjusting the 
optical depth, it becomes much poorer at the other 
end. 

Just as for the synchrotron model, a more sys- 
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tematic treatment of thermal absorption is to solve 
for the spectrum as well as its derivative by finding 
the parameters of best fit. To achieve this, we also 
allow the synchrotron energy spectral index, p, to 
vary freely which, as before, fixes the optically thin 
synchrotron spectral index, a = (1 — p) 12. The 


D2 


amplitude scaling, C, incorporates C of Eqn. 
as well as a term to convert from specific intensity 
to flux density for comparison to observations. Fi¬ 
nally, we combine the remaining unknown param¬ 
eters into a single free parameter, k, which is a 
variant of Eqn. |D1| that also includes the line of 
sight distance, sn (in units of 10 ^^ cm), and the 
frequency at which the synchrotron spectrum be¬ 
comes optically thick, z^i(9), in GHz. That is 




Sl7 




( 10 ) 


where nec(4) is the thermal electron density in 
units of 10 "^ cm“^ and T 4 is the ‘cold’ electron 
temperature in units of 10^ K. 

The resulting best fit parameters are listed in 
Table and are shown as green curves in Fig. 
This fit is demonstrably poorer than that of pure 
SSA, consistent with our statements above. One 
notes particularly the poor fit to the spectral index 
in C band. 


4-2.3. Synchrotron spectrum with a foreground 
thermal screen 

It is natural to imagine thermal absorption 
ocurring in a sheath of mixed cold and hot par¬ 
ticles surrounding the synchrotron jet, given the 
lack of linear polarization that we observe (see Ap¬ 
pendix H- In principle this could be a dynamic 
sheath of cold gas entrained (that is mixed) with 
the jet particles, or it could be be a stationary 
surrounding cloud. 

We continue to assume that the synchrotron 
emitting particles are additionally self-absorbed. 
Then the equations of the fit are 


S^, 

a 


-exp(-^))exp(-^), 

q exp(-A) 


We take q = (2 + p/2), z = v/vi as before, and 
K as above. This model has four free parameters 


p,vi,k and the amplitude C. The results are 
given in table and by the blue curve in Fig. 

The fit is no better than that of SSA alone al¬ 
though it is equally good. The high frequency 
spectral index problem of the pure thermal ab¬ 
sorption has been removed. However this has been 
accomplished by removing almost completely the 
thermal absorption at C band (oc e“ t 6 ), which 
leaves the optically thin synchrotron emission un¬ 
changed there. Nevertheless this is an acceptable 
configuration. 


We observe also from Fig. that any of these 
simple models fails to reproduce the fluxes to ob¬ 
servational accuracy. This is not too surprising 
since we have ignored any real geometrical struc¬ 
ture in the source. For example we see from the 
definition of k that the line of sight required for 
a given absorption is rather sensitive to T4. More¬ 
over a detailed radiative transfer calculation in the 
style of [Jones ' (1988) for a given (partially opti¬ 
cally thick) geometry may be required. 


It is interesting to note that Wang et al. (2014) 


have invoked very similar ideas to our own regard¬ 
ing SSA and SSC (Synchrotron Self-Compton) 
in the description of a more distant TDE (Sw 
J16449.3-1-573451). Their geometry requires two 
nested conical jets of different speeds and total 
energies. The slower outer jet may be compared 
to our idea of a dynamic sheath. However we dif¬ 
fer from these authors in using the jet to produce 
the early X-ray emission in a shock lobe followed 
by the radio emission after the decelerated jet has 
expanded. 


A detailed study is justified for their source 
since much more data in the x-ray and radio is 
available than is the case for NGC 4845. For this 
reason we will adopt our simplest model with the 
smallest number of parameters, namely pure SSA, 
in the next section on variability. Subsequent more 
detailed work after future VLBI obervations may 
be justified. 


4.3. Variability of the Spectrum 

4-3.1. A Single Injection of particles followed by 
fading does not fit the data 


In a paper discussing a distant GPS source ( Ori- 
enti et al. [2010 1 it was concluded that the overall 
source spectrum was best explained by an ensem¬ 
ble of electrons in which acceleration had ceased. 
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The authors called this the ‘fader’ phase. Their 
extensively sampled spectra were fit in terms of an 
‘injection spectral index’ and in terms of the ratio 
of the acceleration cutoff time to the synchrotron 
lifetime at peak frequency. 

In the most successful ht to the overall spec¬ 
trum (that is with a reasonable injection spectral 
index of a « —0.7) they found a cutoff time which 
was about 20% of the synchrotron lifetime of 
several thousand years. However these authors 
did not observe a fresh outburst in this source 
so the expected dependence on time could not be 
checked. This is the distinct advantage of detect¬ 
ing the combined X-ray and radio emission of the 
outburst in NGC 4845. 

We show in Appendixj^that the time evolution 
of the spectra shown in Fig. does not match this 
simple picture. 

4-3.2. Spectral Variability as a Result of Adia¬ 
batic Expansion 

We now wish to understand the time evolu¬ 
tion of the source spectrum, again summarized 
in Fig. (from epochs Ti to T^), which we take 
to be the best representation of the data. As re¬ 
marked above, we assume SSA uniquely since this 
minimizes the number of free parameters and pro¬ 
vides a fit to the observations which is equally as 
good as a synchrotron spectrum with an additional 
thermal screen. We therefore suppose the energy 
spectrum, N{E), of the relativistic particles to be 
evolving due to energy loss through adiabatic ex¬ 
pansion. 

This mechanism presupposes pressure equilib¬ 
rium with the medium surrounding the evolving 
source, against which the relativistic particles are 
working. This medium may be thermal gas in 
a ‘sheath’ around the source, perhaps combined 
with (or indeed dominated by) the ‘hoop stress’ 
of a helical magnetic field. Such a local confining 
magnetic held must be created by an externally 
driven current. 

During adiabatic expansion the equation for the 
total electron density per unit volume and per unit 
energy, n{r,E,t), is (ignoring acceleration, diffu¬ 
sion and particle source terms) 

9*71 + V • (nv) = i(V • v)5£:(nF;). (12) 

where r is a position vector, v is a velocity, and 


t is time. For our purposes we integrate this over 
an expanding volume that contains the same total 
number of particles at any given time, in order 
to get the familiar ‘diffusion-loss’ equation (e.g. 
Longair 
energy interval, N 

dtN{E,t) = ^{W ■v)dE{NE). (13) 

We have assumed that the velocity divergence can 
be expressed wholly as a function of time (includ¬ 
ing a constant) in order to perform the spatial in¬ 
tegration. 

For simplicity (but see Appendix]^ we assume 
that the source expands with the uniform radial 
velocity, VaSr, where the radial direction, r, is 
along the axis of the cone. Hence V • v = 2vofr = 
2lt , since r = Vot. The solution for N{E,t) is 
then found (e.g. by the method of characteristics) 
to be 

/ f \ (2(l-p)/3) 

fV(£;,t) = iVoA-M-j . (14) 

Here, the reference time, ti = 342 days (elapsed 
time since the X-ray peak at to) occurs at the time 
stamp, TI. The energy spectrum at time, ti, has 
clearly been taken as usual to be N{E) = NoE~p. 
Subsequently the spectrum is unchanged but now, 
when p = 2, 

A„(t)=iV„(ti)(t/ti)-2/3 (15) 


for the number of particles per unit 


The frequency at synchrotron optical depth one 
1^1 (which is proportional to the peak frequency, 
Vrn) satisfies (Pacholczyk ||1970 ) 


j/i cx 


(16) 


when p = 2. If we suppose that the line of 
sight size, s cx. t (as does r), and use the form 
No{t) found above, then 

vi oc [t/tiy/^B]^^. (17) 


We would like to explain the ratios of the peak 
frequencies which, from Table[^ are VmiTl)/Vm(T2.) 
4.87/4.03 « 1.2 and Vm{Tl)/= 4.87/3.22 « 
1.5 in terms of this expansion. If B± oc 1/s so that 
the azimuthal field dominates (the azimuthal field 
wraps around the cone/jet. Fig. 1^, then we cal¬ 
culate the ratios as (1 -I- 56/34^/® « 1.1 and 
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(1 + 196/342)^/® « 1.3, respectively. If, at the 
other extreme, we suppose that the radial field 
dominates B±, then B± oc 1/s^ and we cal¬ 
culate the ratios as (1 -I- 56/342)^^/® « 1.2 and 
(1 -I- 196/342)^^/® « 1.7, respectively. Overall, ei¬ 
ther decline in magnetic field approximately fits 
the data; however, we show below that B± oc 1/s® 
is required to explain the amplitude variations. 
This suggests a rather losely wound, conical mag¬ 
netic helix together with moderate particle pitch 
angles in the source. 


Finally, we wish to consider the changing am¬ 
plitudes of these peaks. The specific intensity at 
the peak frequency, is given by (Pacholczyk 
19701 P- 97) 


oc Bj_ 


1/2 5/2 
‘'l 


(18) 


Applying Eqn. nil this becomes 


result in an increasing flux density with time which 
is clearly not observed. Thus, hereafter, we adopt 
Bj_ oc 1/s®. 

For the sake of completeness in the context of 
our model, we can also hnd the time dependence of 
the flux density, for any hxed frequency, v, in 
the optically thin limit (as opposed to since 
I'm changes with time as found above). For this 
case, 

S, = 0,® J, = 0,® Jiz) (21) 


where J{z) is given by Eqn. and is related to 
tiy a simple constant. T he time dependence 
of 6s^ is given by Eqn. and, in the optically 
thin limit with p = 2, J{z) oc Since vi oc 


t ^^/® as already found from Eqn 17 with B± 
1 /s® and s oc t, we find finally 


oc (22) 


oc (iVo(t))®/®B_L^/®s®/® oc (19) 


The latter time dependence results from the fact 
that each of the above terms can be expressed as 
a function of time using Eqn. 15 B± oc 1/s® from 


above, and taking a constant velocity so that s oc t. 


It remains to convert from specific intensity to 
the observed flux density which requires consid¬ 
ering how the projected area of the source varies 
with time. For spherical outflow, this dependence 
would be oc s® oc t®. However, a jet/cone 
geometry should have Vr >> v±, where v± 
are velocities along the axis of the cone (the radial 
direction) and perpendicular to it, respectively, in 
which case 0^^ oc s oc /. In the unlikely event 
that t = 0 (the line of sight is directly along the 
axis. Fig. [^, then 0^® « constant for a narrow 
outflow. Adopting the jet/cone geometry, then 
the time dependence for the peak flux density is 
finally 


= 0«®/^^ oc t-i®/i8 « l/t (20) 


Again, from Table the amplitudes of the 
peaks at time stamps, T1 to T3, are in the ratio, 
St 2 /Sti = 398.5/432.7 = 0.92, and Sts/Sti = 
275.4/432.7 = 0.64, and the predicted ratios are, 
(1-f 56/342)-! ^ Q gg (1-H196/342)-! = 0.64, 
i.e. in agreement to within 7% or better. 

If we repeat the above, assuming that B± oc 
1 /s, then the result is oc t^/® which would 


Such a dependence might be expected at C- 
band where the emission is optically thin. Note, 
however, that we base our analysis on the (exact) 
polynomial spectral fits to all of the data (fluxes 
and slopes) in Fig. rather than on partial spec¬ 
tral fits to the data as displayed in Fig. |^for T2. 
It is readily seen that the shape of the polyno¬ 
mial curves is not always consistent with the C- 
band flux density lying in the optically thin re¬ 
gion, where the decline should apply. This 

is because the value oi z = v^jv^ in C-band for 
T1 is found to be only 1.74. Fig. then shows 
that a% = —0.224. For T2, z = 2.1 and thus 
ttQ = —0.341. That is, with increasing time, the 
spectrum becomes more optically thin, consistent 
with what might be expected from adiabatic ex¬ 
pansion, but the fully optically thin regime is later 
than T2. 

Although the T2 polynomial spectrum is still 
not convincingly optically thin, it is the best point 
from which to try the extrapolation that was 
carried out for the C-band curve in Fig. (dashed 
curve) as an alternative to the simpler linear ex¬ 
trapolation (dotted curve). 

In summary, both the changing frequencies and 
amplitudes of the radio spectra are in reasonable 
agreement with a self-absorbed cone/jet model 
which is adiabatically expanding and originates at 
the time of the hard X-ray flare peak. 
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4.4. Circular and Linear Polarization 

The most evidently remarkable results concern¬ 
ing polarization are reported in Tables and 
NGC 4845 is in the extension of the Virgo clus¬ 
ter at high Galactic latitude (6 = -1-74°) so these 
measurements are unlikely to be contaminated by 
the Milky Way. 

The first table reports significant values (be¬ 
tween 2 and 3%) of circular polarization (circularly 
polarized fraction, or cpf) in L-band at all reso¬ 
lutions. These are large compared to maximum 
values previously observed in GPS sources (« 1% 
e.g. O’Sullivan et al. 20131, although not highly 
discordant, given the uncertainties (Sect. 3.21. 


Equally remarkable is the non-detection of circu¬ 
lar polarization at C-band. There is evidently an 
extremely rapid variation of circular polarization 
with frequency, reported in the second table, i.e. 
the in-band spectral index of the circularly po¬ 
larized flux lies between —2.2 and —3.4. This be¬ 
haviour is coupled with no reliable detection of lin¬ 
ear polarization. We need to constrain the source 
structure to fit these facts. 

A satisfactory explanation of our observed po¬ 
larization behaviour is difficult to find in the lit¬ 
erature. In this section, we investigate several 
approaches, and then discuss our proposed expla¬ 
nation which involves conversion from linear po¬ 
larization to circular polarization via generalized 
Faraday rotation (see also O’Sullivan et al. 

The intrinsic circular polarization percentage of 
optically thin synchrotron radiation is cot(()/ 7 e, 
where (f) is the pitch angle and 7 e is the elec¬ 
tron Lorentz factor corresponding to the peak fre¬ 
quency. For a 0.01 G field, 7 e « 200 in L-band 
(Eqn. |F1| , or 7 e « 100 for a 0.04 G field. Then 
one would need </> = 14° —>• 26° (for 7 e « 200 —>■ 
100 , respectively) in the observer’s frame to ac¬ 
count for a 2% circular polarization percentage. 
Such a pitch angle reduces the emissivity to about 
12 —> 30 % of the maximum (cx {B 
which, to maintain the observed flux at the same 
frequency, would require a rather large magnetic 
field (10“^ G) and hence magnetic energy (^ 10"^® 
ergs). 

The major difficulty with this explanation, 
however, is that the intrinsic polarization predicts 
a spectral index (ay) for the circularly polarized 
flux density of —0.5. This conflicts with what we 


(Table an d others (see discussion in O’Sullivan 
et al. |[20l3 ) have observed. A second difficulty is 
that we detect little or no linear polarization. This 
would not be expected for optically thin, homoge¬ 
neous emission. In any case we have already made 
a case that our frequency observations embrace 
the transition to self-absorption (Sect. [4(2T ). The 
optically thick regime is also necessary for our 
proposed explanation by generalized Faraday ro¬ 
tation (e.g. Pacholczyk [1973 


O’Sullivan et al. 20131. 


Jones 1988 Ken- 


nett & Melrose 1998 [Zheleznyakov fc Koryagin 


2002 Macquart]2002f Beckert fc Falcke 2002 


Jones (19881 gives detailed polarization calcu¬ 


lations for a generalization of a definite jet model 
(see Blandford fc Kbnigl ||1979) and finds circular 


polarization levels in the simulations comparable 
to our own. Unfortunately the calculated linear 
polarization tends to be of the same order, and 
Jones does not discuss the frequency dependence 
in the optical depth transition zone explicitly. His 
analytic discussion focuses on the surface where 
the optical depth is one, whereas we suspect the 
essential behaviour to be in transition across this 
value. The results are complicated by the assump¬ 
tion of inhomogeneity comprised of turbulent cells. 


eke ( 

20021 . 

In 

Pacholczyk 

(1973 


circularly polarized flux in a source comprised of 
mainly cold electrons with an admixture of rela¬ 
tivistic particles is discussed. However ay = — 1 in 
this limit. Linear polarization is suppressed only 
by the average effect of the turbulent cells, if the 
field is indeed randomly distributed and the Fara¬ 
day rotation per cell is large. When the Faraday 
rotation per cell is small, the circular polarization 
has the dependence near optical depth unity. 


but the linear polarization is larger (Beckert & Fal¬ 
cke 1200 ^ . This may have to be depolarized in an 
inhomogeneous, Faraday rotating screen. 

4-.4-1- Generalized Faraday Rotation 


Generalized Faraday rotation occurs (e.g. Ken- 


nett & Melrose 19981 when the orthogonal nor¬ 


mal modes in the source plasma are elliptically 
polarized rather than circularly polarized (effec¬ 
tively linearly polarized) as is the case for stan¬ 
dard Faraday rotation. The latter process occurs 
in a cold (thermal) electron plasma and is pro- 
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portional to -By. The former process requires 
an admixture (becoming dominant in a relativis¬ 
tic limit) of relativistic electrons and is propor¬ 
tional to B\. Normally the generalized process 
produces a cyclic generation of circularly polar¬ 
ized flux from linearly polarized flux (this is the 
‘conversion’). That is the case, provided that the 
plane of linear polarization rotates along the line 
of sight so that there is appropriate rotation on 


the Poincare sphere (Kennett & Melrose 19981. 
Thus, either the magnetic field rotates or some 
cold plasma must be present to Faraday rotate the 
plane of polarization. The cyclic variation modu¬ 
lates the frequency dependence of the cumulative 
phase shift between the two modes. If the wave 
number difference is Afc then this phase shift is 
Afcs, where s is the line of sight distance through 
the source. 

It is important to note that the steep in-band 
V spectral indices of Table apply over a rela¬ 
tively narrow frequency range (500 MHz) where 
the transition to self-absorption is occurring, as is 
indicated by the spectral index awl. In a region 


where relativistic particles dominate (Kennett & 


Melrose 19981 and optical depth effects are ne¬ 


glected, the percentage cpf, V/S, is proportional 
to Afc oc which is roughly consistent with our 
observations. As the relativistic particles become 
less important relative to the cold electrons (a cri¬ 
terion is given in Kennett & Melrose 19981 the 


variation tends to i/ ^ ( jPacholczyk [1973 1. 

We can consider the dominance criterion for 
relativistic particles (denoted with subscript, er) 
over thermal (‘cold’) particles, (denoted, ec). This 
criterion is n^c 7 ifT-er In (27rz//z/gB7o) (recall 
Vg is the gyrofrequency, Eqn. [^, which requires 
Uec ^ 2.3 X 10^ cm~^ if B = 0.01, 70 « 5 (e.g. 
Beckert fc Falcke~||2002 1 at L-band. This is about 


the same order of magnitude as the density of rel¬ 
ativistic particles estimated in Sect. [T2T] and im¬ 
plies that the thermal mass in the source is prob¬ 
ably of the same order of magnitude as the rela¬ 
tivistic particle mass. We will pursue this further 
in Appendix [E} 

The only similar result of which we are aware 
is for the quasar PKS B2126-158 as reported in 
O’Sullivan et al. (2013). They find ay to be 


—3 ± 0.4 at somewhat higher frequency (8 GHz) 
compared to our result. Their brightness spec¬ 
trum peaks at a frequency similar to what we 


infer in Fig. Moreover at lower frequencies, 
where oc zz as in our spectrum, they detect 
rising circular polarization between 0.5% and 1% 
in what is the optically thick region. This quasar 
source thus displays a slightly broader transition 
between synchrotron self-absorption and an opti¬ 
cally thin regime. Significantly when compared to 
our results, the observed linear polarization is also 
nearly zero near the peak of the spectrum. This 
is where their observed circular polarization is in 
fact largest, much as in our case. These authors 
explain their —3±0.4 spectral index in circular po¬ 
larization in the optically thin regime, via linear 


to circular conversion (Kennett & Melrose 1998) 


(see next subsection) which we also adopt in part, 
although our measurements are taken closer to the 
optically thick transition. 


4-4-^- Conversion of Linear to Circular Polar¬ 
ization 


The mechanism described above, of course, re¬ 
quires an initial linear polarization to be rotated 
on the Poincare sphere. Fortunately, [Beckert 


Falcke (2002) give a relatively simple set of equa¬ 


tions together with a complete set of absorption 
and emission coefficients for a relativistic plasma, 
which allow us to investigate the spectral depen¬ 
dence and conversion from linear to circular po¬ 
larization in more detail. We develop these in Ap¬ 
pendix and extend their analysis to include the 
case p = 2 , as implied by our total intensity spec¬ 
trum (Sect. [T2T ). 


We consider the two limiting cases of small and 
large Faraday depth. For the core of NGC 4845, 
we also show numerically that such conversion can 
explain the circularly polarized flux that is ob¬ 
served. However, the true condition in the source 
is likely to be intermediate between these limits. 
The limit that is closer to what we observe is cer¬ 
tainly the low Faraday rotation source plasma. 
The circularly polarized flux is given by Eqn. |E11[ 
as a blend of synchrotron emissivity and rotational 
conversion, which agrees with our observed ay. 
Should B(—2) be as large as 10, the required linear 
polarization for conversion is reduced to typical 
values (3 or 4 %). To reduce the linear polariza¬ 
tion further, however (i.e. to less than 0.5% as 
observed), it must still be suppressed by an inho¬ 
mogeneous, depolarizing medium. The presence of 
such a medium is not unusual for edge-on galaxies; 
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a comparison to other galaxies in the CHANG-ES 


sample (Wiegert et al. [2015) shows that a large 


number also lack linearly polarized flux at L-band. 

We conclude that to explain our observations 
we need only that the relativistic particles are 
present and are dominant at L-band. The bound¬ 
ary conditions U = Uo and Q = E = 0 do suggest 
that we are observing an organized synchrotron 
source through an inhomogeneous, medium of 


& Koryagin 

2002) 

source near t 

re jet 


For a coherent synchrotron 


general scheme), then we may suppose that our 
calculations apply to a sheath where relativistic 
particles still dominate. The depolarizing screen 
around the jet possessing dominant Faraday ro¬ 
tation may gradually appear due to an influx of 
entrained thermal particles. The circular polariza¬ 
tion incident on the screen should not be affected 
by the Faraday rotation in the screen. This pic¬ 


ture is similar to a two jet model (Wang et al. 


20141, except that the outer jet here contains a 


mixture of relativistic and thermal electrons. The 
thermal particles presumably become dominant at 
larger scales perhaps due to entrainment of a sur¬ 
rounding thermal cloud. 


4.5. A Proposed Radio/X-ray Connection 

4-. 5.1. Inverse Compton Radiation 


We first consider whether inverse Compton ra¬ 
diation (the so-called synchrotron self-Compton, 
or SSC mechanism, for up-scattering from radio 
to X-ray energies) can explain the X-ray emis¬ 
sion, given the properties as measured in the radio 
source now. 

Inverse Compton emission is related to syn¬ 
chrotron emission from an ensemble of rela¬ 
tivistic particles through the ratio of the radia¬ 
tion energy density to the magnetic energy den¬ 
sity, Urad/Umag, that is, the SSC flux would be 
Sssc = Sn Urad/Umag, where Sr (ergs“i cm“2) 
is the measured radio flux. We estimate the ra¬ 
diation energy density from Urad = (Att/c)/, so 
that 


Urad 


16 vSa 


5x10-"%®^ 

()s\j^o.sy 


ergs cm 
(23) 


We have taken the maximum from the earliest 
epoch in Fig. when the frequency is 4.87 GHz 


(letting ;/( 4 . 87 ) = 1) and the flux density at that 
frequency is 433 mJy (Table [^, or 5(4,87) = 1- 
Using Umag = -B^/(87r), we Arid the ratio 


Urad ^ >S'(4.87)i^(4.87) 

Umag ' 9sijnasYB{-2Y ' 


(24) 


We might therefore expect inverse Compton emis¬ 
sion during the radio outburst at about 10 per¬ 
cent of the radio flux density using B{—2) = 1 
(0.01 Gauss) for the field strength (Appendix]^. 
This emission should be in a band about the peak 
value ')/vpeak. Using our previous estimate for 7 e 
(estimated for L-band, Sect. 4.41 we expect this to 
lie in the near infrared to optical bands - not the 
X-ray. 

To determine the expected X-ray flux, we can 
find the maximum radio flux from our observations 
by integrating over the band from 0.5 to 10 GHz 
for the highest curve of Fig.|^ finding « 2 x 10“^^ 
erg s“^ cm“^. The X-ray light curve in the 17.3-80 
keV band ranges from 10“^^ —>■ 8 x 10“^^ 


erg s 


cm 


( Nikolajuk &: Walter ||2013 their Fig. 8). 
The last point, in particular, in Fig. 8 of |Niko- 


lajuk & Walter (2013) overlaps in time with our 


observations. Whatever the origin of this X-ray 
emission is, it cannot be due to SSC from the 
radio. Not only is the luminosity low, but the 
Compton scattered frequency of the radio photons 
falls far below the X-ray band as indicated above. 
Thus there is another source of X-ray emission in¬ 
dependent of the SSC mechanism. In |Nikolajuk| 
& Walter (2013) it is recognized that this last 


point does not fit the expected tidal disruption 
behaviour that is found near the peak X-ray emis¬ 
sion. It is likely that some late behaviour arising 
near the accreting disc of the system (which is per¬ 
haps at ^ 1 AU from the black hole) is responsible. 


4 . 5 . 2 . Inverse Compton Radiation from an Evolv¬ 
ing Source 


If the radio emission does not explain the X-ray 
emission via SSC for current conditions, can we 
explain the X-ray emission by extrapolating the 
radio conditions to the epoch of the peak X-ray 
emission? 


A remarkably similar set of observations to our 
own were reported for the radio emission associ¬ 


ated with a Swift Gamma-ray source (Zauderer et 
2011|. These authors concluded that they 


al. 
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were seeing the birth of a parsec scale radio jet 
in the core of a rather distant galaxy. Like us 
(previous subsection), these authors rejected in¬ 
verse Compton emission because of the luminosity 
shortfall, and also because of the lack of correla¬ 
tion between the radio and X-ray variability. 


However, in our case, the X-ray energy distribu¬ 
tion power law near the peak emission (px ~ 2.2, 


Nikolajuk &: Walter 2013) is similar to the radio 
energy distribution power law (« 2, Sect 4.2.1) 
and both X-ray and radio emission consist of a 
single strong point-like source at the center of the 
galaxy. We are thus further motivated to attempt 
to establish a potential link in the jet development. 
A complete description of our development is given 
in Appendix [F| and is within the context of the jet 
model of Appendix We summarize the argu¬ 
ments and results here. 


Recall that timescale arguments (Sect. [4d| ) sug¬ 
gest a radio source size of order 0.1 pc (~ 1 mas) 
and, if the jet is relativistic (Appendix]^, then 
we expect that the shocked X-ray emitting region 
at the earlier time (the epoch of the peak X-ray 
emission about 1 year earlier) is an order of mag¬ 
nitude smaller in size, a region that we refer to as 
the X-ray ‘lobe’. The observed flux density of the 
circular polarization is consistent with a magnetic 
field strength of order 10“^ G (these are order of 
magnitude estimates). Along with the observed 
flux density of the peak of the radio emission at 
time, Tl, these are the only inputs to our model, 
the inputs being constrained by observation. 


The extrapolation allows the magnetic field 
strength in the X-ray lobe to be estimated since 
the variation with size is known for a constant ve¬ 
locity jet [B oc 1/r^, Appendix [A|) . Since the X- 
ray frequency is known, an estimate of the Lorentz 
factor that is required to produce the X-rays via 
SSC can be determined along with the critical 
frequency of synchrotron photons in the X-ray 
source. We find that the synchrotron emission in 
the X-ray lobe would peak in the infra-red. 


The flux of the IR synchrotron radiation in the 
X-ray lobe can be found by relating the currently 
observed synchrotron emission of the radio source 
to the IR synchrotron emission in the X-ray lobe, 
knowing how synchrotron emission depends on the 
Lorentz factor, the magnetic field, the frequency, 
the source size, the electron energy power law 
(p = 2) and the electron particle constant, Nq, 


(the latter quantity depends on time in a known 
fashion for an adiabatically expanding source, see 
Sect. Its I. The X-ray flux from SSC, 


S 


_ XssC! 

the X-ray lobe can then be determined from the 
IR synchrotron emission and magnetic field (pre¬ 
vious subsection). 

Our calculations (Appendix for the above- 
indicated conditions, show that Uradm/Umag = 
1.1 and Sxssc — ^-O x 10“^^ erg s“^ cm“^, the 
latter value falling short of the peak obser ved X- 
ray flux of Sx = 8.0 x 10“^^ erg s“^ cm“^ (Niko¬ 


lajuk & Walter 2013) by a factor of 4. The re¬ 


sult is, however, very sensitive to the magnetic 
field in the X-ray lobe. We have, so far, allowed 
for an increase in the magnetic field to the ear¬ 
lier epoch that only accounts for a change in size 
via jet geometry, and have not yet allowed for the 
possibility that shocks in the X-ray lobe may en¬ 
hance the magnetic field further. If we include a 
modest enhancement of the magnetic field (a fac¬ 
tor of 1.8, whereas a strong shock would give 4), 
then we find the same ratio of Uradm/Umag and 
Sxssc = 8 X 10“^^ erg cm“^, in agreement 
with the X-ray observations. 

Given these results, we are encouraged to dis¬ 
cuss such a SSC X-ray source in more detail (see 
Fig. [T]). We imagine that largely thermal mate¬ 
rial at a very high temperature is expelled from 
near the black hole due to the tidal disruption of a 
Jupiter-sized object. The disrupted material also 
forms a disc or torus near 1 AU. The ejected ther¬ 
mal material is collimated into a fast beam by this 
disc. Subsequently the beam collides with a cloud 
of interstellar material at about 0.01 pc and forms 
an X-ray-emitting region (the X-ray ‘lobe’). This 
lobe consists of a forward shock in the cloud, a re¬ 
verse shock in the beam, and a separating contact 
discontinuity in the medium between the shocks. 
The model is much as was described in another 


context in Henriksen, Mirabel & Ptuskin (1991). 


For the shock wave acceleration to be effec¬ 
tive in accelerating the electrons, the acceleration 
time, tacc, must be less than the synchrotron life¬ 


time. In Henriksen, Mirabel & Ptuskin (1991) 


this was estimated in their equation (7). For rela¬ 
tivistic shock speeds (it is sufficient for vshock to be 
> 0.001c) and our magnetic held {B in the X-ray 
lobe of 1.8 G), then tacc < 10 min (proportion¬ 
ately shorter for higher vshock)- The synchrotron 
lifetime of an electron (Eqn. Cl) under these con- 
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ditions with 7 e = 1109 (Appendix[F| on the other 
hand, is about 20 minutes. Hence the acceleration 
time is sufficiently short. 


This result also indicates that there must be 
continuous injection or acceleration of the elec¬ 
trons in the X-ray lobe over a time during which 
the SSC mechanism is producing the hard X-ray 


emission. From the X-ray light curve (Nikolajuk & 


Walter 20131, this would occur around the peak of 
the light curve and up to « 120 days after the peak 
if the first cluster of X-ray measurements after the 
peak is produced in the same way. After that 
time, the X-ray light curve declines precipitously 
(suggesting that there is no longer continuous in¬ 
jection/acceleration) revealing a secondary X-ray- 
emitting source, possibly associated with the ac¬ 
cretion disk. 


If the reverse shock is a very strong shock, it 
produces an accelerated distribution of electrons 


with energy power law, NoE~'^/'^ (Maikov & Drury 
2001), i.e. the injection spectrum is described 
by Pinj = 3/2. If this is the continuous injec¬ 
tion spectrum for the electrons from this shock, 
then according to Longair (19941, above an en¬ 
ergy given by 


E., = 


1.2 X 10^ 


BlU(l 


-eF 


(25) 


{ts is the source lifetime in years and is the 
field in Gauss), the spectrum steepens from syn¬ 
chrotron losses such that p = pmj -1-1 = 2.5. 
The required electrons are well above this energy 
and therefore we expect a power law, E~^'^, from 
this shock. If the forward shock is less strong, 
then we expect the test particle result, p = 2, to 
be given to the electrons by this shock. The elec¬ 
trons scattered upstream of the forward shock will 
form the radio spectrum with p = 2. These must 
be a significant fraction of the electrons scatter¬ 
ing photons to X-ray energies, but not in energy 
space. These accelerated electrons that leave the 
shock in the forward direction cool rapidly by syn¬ 
chrotron radiation and eventually form the adia- 
batically expanding radio lobe that we infer from 
our observations. The cooling time is rapid enough 
to convert the high energy electrons that leave the 
X-ray lobe into the low energy radio electrons that 
are observed about a year after the peak. 

In summary, in the context of this model, the 
regions of the shocks, of size « 1000 AU and mag¬ 


netic field « 1 G can produce X-ray emission via 
SSC (the X-ray ’lobes’) with outflow in the for¬ 
ward direction developing into the radio jet. Rel¬ 
ativistic electrons cool rapidly after leaving the 
shocked region, and develop into an adiabatically 
expanding outflow. 


5. Conclusions 


Our main conclusions are as follows. 

1. The radio emission from the Virgo Cluster 
spiral, NGC 4845, is dominated by a single unre¬ 
solved point source that has increased by a factor 
of « 6 since the 1995 NVSS survey and also varied 
over the « 6 month time scale of the CHANG-ES 
observations (Table |^. Prior to NVSS, there is 
no evidence for radio variability from the few data 
points that exist (Sect. [3X4 ). However, some ver¬ 
tical radio extensions as well as a distinct optical 
cone (Fig. suggest that previous outflow has 
likely occurred. 

2. A small central disk of diameter, 1.8 kpc, 
surrounds the nucleus but constitutes only 19 mJy 
at L-band (less than 10% of the total flux density) 
and 7 mJy at C-band (less than 2% of the total). 
The spectral index of this disk is adisk = —0.74 
(Sect. |3n| ). 

3. The CHANG-ES observations overlap in 
time with the hard X-ray light curve of [Nikolajuk] 
& Walter (2013 their Fig. 8). The X-ray flare, in¬ 
terpreted as being due to a tidal disruption event 
(TDE), provides a well-defined reference point for 
the origin of the current radio emission. We pro¬ 
vide a simple jet model for the radio emission (Ap- 
pendix[A|and Fig. [ 7 ]) which is relevant whether the 
X-ray outflow originated with a TDE or is due to 
infall from an accretion disk, though the former 
may lead to the latter in any case. 

4. Variability suggests a source size of 9s < 
4sin (t) mas (of order 0.1 pc), where l is the view¬ 


ing angle (Sect. 4.1). 

5. The wide bands used in the CHANG-ES sur¬ 
vey have allowed us to determine in-band spectral 
indices and to solve explicitly for the spectrum of 
the AGN as a function of time (Sect. [3.1.5 and 
Fig.|§. At L-band, ai.sy qhz ~ -1-1.0 and at C- 
band, ae.oo GHz ~ —0.5 (Table |^, suggesting an 
energy spectral index for the relativistic electrons 
of p = 2.0 and that NGC 4845 is a nearby exam¬ 
ple of a ‘GigaHertz-peaked spectrum’ (GPS) radio 
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source. 

6. The turn-over in the spectrum can be ex¬ 
plained by synchrotron self-absorption (Sect. 4.2.1). 
However a contribution from thermal absorption 
at L band is possible. 

7. The peak of the radio spectrum both declines 
and moves to lower frequencies with time (Fig.|^. 
This behaviour can be understood in terms of adi¬ 
abatic expansion (Sect. 4.3). 

8. We detect circular polarization (Sect. [4(4| ) of 
order a few % in the source, although no signifi¬ 
cant linear polarization, at L-band (Table [^. The 
circularly polarized flux has a steep spectral index 
(» —3, consistent with other authors) and can be 
explained by conversion from linear polarization 
- a natural consequence of generalized Faraday 
rotation for a source that is dominated by rela¬ 
tivistic particles (see Beckert fc Falcke~||2002 and 
Appendix . A further reduction of linearly po¬ 
larized flux requires an additional foreground de¬ 
polarizing screen. 

9. Inverse Compton emission from the currently 
observed radio source is insufficient to account for 
the X-ray emission detected by|Nikolajuk & Wal-| 
ter|(|201^. 


10. However, if we extrapolate to the peak of 
the X-ray light curve at which time the source size 
is smaller and the magnetic held higher (Fig. [^, 
the X-ray emission can be explained by inverse 
Compton upscattering by the relativistic particles 
(SSC), a result supported on energetic grounds 
and also by the similarity between spectral indices 
at radio and X-ray wavelengths (Sect. [T^). 


NGC 4845 appears to be a young nuclear radio 
jet/outflow and is, by far, the closest known GPS 
source. As such, it provides a unique opportunity 
to study an evolving radio source associated with 
a measured hard X-ray burst. The radio source, 
after 5 years at an assumed bulk outflow speed 
of c/-\/3 (Appendix [a|, would have expanded to 
about 1 pc, or 12 mas. This size is within the 
realm of Very Long Baseline Interferometry. With 
observations over several years, it may be possible 
to directly observe the proper motion and evolu¬ 
tion of this source, and further piece together its 
connection to hard X-ray outflows. 
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A. The Standard Jet Model 


We are imagining a shocked section of a conical outflow (or possibly the pseudo core discovered in Jones 
19881 as the source region. Fig. shows the geometry. 

The jet model summarized here has been assumed as the synchrotron source in AGNs by many authors, 


beginning with Blandford & Konigl (1979) and then Jones (1988). However the dynamical details of such 


sources are not readily found in the literature. Moreover we are able to suggest a novel observational idea 
based on the net current flow through the jet. For these reasons we have decided to include a description 
of the model in this appendix. Ultimately it would form the basis of an analysis in the style of |Jones~| 
(1988), which is the level required to match more abundant observations. The simple model described here 


is similar to that used by Jones (1988). However we treat the magnetic field in more detail and we include 


a decelerating branch of the outflow. Whether this branch or the more familiar accelerating branch appears, 
depends on ambient boundary conditions. In addition, we consider the possibility of change in the direction 
of the current flow in the jet and corresponding change in the magnetic field. 

One assumes that the jet originates near the compact object as a gas dominated by relativistically hot 
particles. Moreover if p is the gas pressure and p is the total energy density (including the thermal energy 
then we assume for simplicity the ‘polytropic’ relation (we take units with c = 1, e.g. Og c is the sound speed 
in usual units) 

P = alp. (Al) 


One defines the conserved rest mass density pt according to p + p = pi{dp/dp), so that after inserting the 
pressure and integrating 

(A2) 


P _ ^ P- y+at 

Ps ps 


The constants with subscripts s are evaluated at some convenient reference point. 

It is also useful, when expressing the dynamical equations, to define the internal energy per unit rest mass 
energy as 

n n ^ n. 11 o 

(A3) 


P Ps Ps 


Hence also 


P = 


jp 

1 + a? 


(A4) 


We assume strictly radial flow, u, in the jet together with a steady state. This renders the transverse 
equilibrium equation useless as it only defines some unknown confining force in the theta direction. We are 
left with energy conservation along the flow direction together with a rest energy conservation law . These 
yield respectively, after a slight rearrangement, when written in the observer frame of reference 


^7('u) = £, (A5) 

and 

r^pj{u)vr = A. (A6) 

Here £ is essentially the constant specific energy in the jet and A is the constant rest ‘mass’ flux through 
the jet. 

For the energy conservation to take this simple form, there can be no Ohmic dissipation or other working 
on the jet by the electro-magnetic field. Assuming that the electric field vanishes in the co-moving frame 
(E = —u A B) and that the magnetic field in the observer frame is B = then E = uB^eg and 

so we require the poloidal current density jg = 0 in order that the dissipation j • E = 0. In fact one expects 


®We retain ‘standard’ nomenclature that is consistent within this Appendix but may differ from the main text, e.g. p refers to 
pressure here whereas it was used as the energy spectral index in the text (likewise for subsequent appendices). 
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to have j || B, although this is not consistent with purely radially moving charges. It can only be the bulk 
flow that moves wholly radially. 

These last two equations in this simplified form suffice to determine the jet dynamics as a function of r, 
as we will report below. However the question as to the evolution of the magnetic field arises. 

The rest mass conservation law and Faraday’s law can be combined to give the equation satisfied by the 
convected magnetic field on each stream line in the form 


-f —1 = 

dt 

which for purely radial velocity integrates to give 

Bj. 

Bfj, _ 

M7(m) 


(m7(u) “ 

/ B^\ u{r) 

Us 



(A7) 


(A8) 


Magnetic flux conservation is contained in these equations provided that V • = 0. 

To understand the significance of B^s we consider the radial component of the current. Restoring c for 
the moment and taking account of the axial symmetry one finds at any radiuf^ 


2 A(6») 

cr sin0 ’ 


(A9) 


where /^(d) is the total current inside the angle 9. For a zero net current through the central source, counting 
both the jet (polar angle 9j) and a likely sheath (polar angle 9sh > 9j), we must have Ir{9sh) — 0. The only 
way in which the azimuthal magnetic field can reverse sign however is if there is a net current towards or 
away from the central source. 

For example let us assume that the radial current density through the jet is of the form 

jr = ^ (tan — tan nO) sin 0, (AlO) 


where tand_|_ is a convenient parameter and n is an integer that determines how narrow is the jet. Then the 
radial current is found to be (n = 2) 


Ir 


2ttI 



sin 20, In (cos 20) 1 

— ) + ^^ + j(l-cos29) 


(All) 


When 0+ = 7r/6 and n = 2 the current passes through zero at 0 = 19.7°. Beyond this angle the net current 
is negative. A balancing current must flow in the more diffuse surroundings or on the surface in order to 
prevent the accumulation of central charge. Such shearing magnetic fields (e.g. A91 can lead to circular 
polarization (e.g. [Zheleznyakov fc Koryagin~ 20021 and possible rapid variations in the observed rotation 
measure (George Heald-private communication). 


The interaction of the radial current density and the azimuthal magnetic field leads to a ‘pinching’ or 
collimating force. The parameters might be chosen so as to balance the transverse equation at the reference 
point or ‘nozzle’. However the assumed conical geometry is not generally consistent with the subsequent 
transverse force equation, so we do not attempt this here. 


The model does not yet allow for the acceleration of the jet by the magnetic field. The Poynting flux in 
the radial direction due to the wound-up ‘spring’ component of the field at any point Vg is u{rs)B‘^^/A'K. 


‘Radius’ refers to the radial direction; the cross-sectional radius of the cone is referred to as r sin d. 
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The associated energy density is then . One might try to include this effect by adding to al in 

the expression for pressure, the quantity However we shall find that this quantity is 

not constant in radius, so the approximation can only apply over a small section of the jet compared to its 
length. We will also ignore this possibility in what follows, but it can be restored locally simply by modifying 
al to al + al- 

The jet velocity profile follows from Eqn. |A5| and Eqn. |A 6 [ once p, is eliminated in favour of ^ and hence 
7 . We can eliminate most of the parameters by writing x = t/R, where (after restoring c) 


PsC V fJ-sS ) 


(A12) 


The velocity in the jet as a function of scaled radius in units of c is then found from 




1 


/3(l-/32)- 


(A13) 


This expression gives the velocity at x for any particular value of although it contains some surprises. 
As X ^ 00 one can have two branches for < 1. In the accelerated branch /3 — >■ 1, while in the decelerated 
branch /3 —^ 0. The choice is made by the boundary conditions at infinity, which apply a finite pressure in 
the decelerated jet flow and vanishing pressure in the accelerated jet flow. 

It is also clear that there is a minimum radius at which these two branches can exist. If we calculate 
dx'^/dp = 0 we find that at this minmum, which we may take to be Xg, 


P 7 



0^(1 - a 2) 


(A14) 


The minimum point represents therefore the origin of the jet in a ‘sonic’ launch point. The two branches 
meet at this point and there is no flow extending to smaller radius with /3 > 0. This behaviour is due to the 
absence of a restraining force such as gravity. The presence of gravity would turn this minimum point into 
a true choke point in the jet flow or wind. Once Eqn. |A13| is solved we find p from Eqn. |A5| and hence /i, p, 
Br and B^ from the various relations. 

We can be more explicit in the interesting case where ap. = 1/3 (that is, the jet vel ocity is the sound 


speed expected for relativistic internal energy, u = c/\/Z, when c is restored), 
parameter-free cubic equation 

/ 3'-/3 + 4 = 0 


Eqn. 


A13 


becomes the 


(A15) 


and Xg = 7^/2 and Ps = l/-\/3. The discriminant of this equation is A = —4 -f 27/x"‘ = 4(—1 -|- x^/x"^). 
When A > 0 so that x < Xg there is only one real root and it is negative. At x = Xg we have A = 0 and 
there is only one real (double) root Ps = 1/ v^. 

In the region of interest, x > Xg, we have A < 0 and there are three real roots but only two are positive. 
The positive roots take a convenient form; for the accelerated branch as 


P 


acc 




arccos 



(AI6) 


and for the decelerated branch as 


2 Un 
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arccos (— 



(A17) 
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These give respectively Pace —>■ 1 and Pdec —>-0 as a;—>- 00 . Ata; = 
Ps = l/v^ as expected. The general behaviour is illustrated in Fig. 
at a;,. 


10 


Xs each expression yields the root 
The curves actually meet vertically 



Fig. 10.— The upper curve is the accelerated jet universal velocity behaviour when a1 = 1/3 in terms of 
X = r/R. The lower curve is the decelerated jet velocity under the same conditions. The pressure at infinity 
goes to zero to establish the upper curve, while it retains a finite value on the lower curve, the curves meet 
vertically at the central value of u = l/\/3 = 0.57735. 

The behaviour at large x is easy to determine using Eqn. |A16| and Eqn. |A17[ For the accelerated jet one 
finds in lowest order (x means equals in the asymptotic limit) 


P 

1 
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1 - 
x, 
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(A18) 


Eqns. A8 then show that asymptotically Br oc x~'^ and oc x~^ as is normally assumed. The behaviour 
of the magnetic field near the source is however more complicated. This is the branch describing the jet 
expanding into a low pressure ambient medium. 

For the decelerated jet the behaviour follows from Eqn. |A17|as 


P 

1 


y/\ - l/a;4’ 


(A19) 
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and the expressions for fj, and p follow from the dependences on ^ given in Egns fAl^ . Boundary conditions 
at infinity are determined by condtions at the origin Xs and vice versa. The radial field varies again as x~'^ in 
lowest order, but the azimuthal field increases as (x {x — \/x^). Even in the absence of shocks, this can 
lead to an outer brightening of the jet. However a dense medium at infinity such that ^oo = £ is required 
to decelerate the jet. This dense medium may be comprised of material swept-up from the ambient medium. 
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B. Source Proper Frame Properties 


We have used quantities in the local rest frame for our calculations. This ignores the possible relativistic 
motion of the source. In the presence of rest frame gas against which the source is expanding, a deceleration 
to non-relativistic source velocity is possible (see e.g. Appendix]^. However it is important to be aware of 
the modifications required by relativistic source motion. We discuss these in this section. Since NGC 4845 
is not at large cosmological redshift, we do not include this effect. 


In the first instance we have argued generally in terms of the standard isotropic, energy power-law, 
distribution of source electrons (or possibly positrons). However a particle moving in the source frame at an 
angle 6' to the relative velocity, appears at an angle 9 to the velocity in the local rest frame where 9 is (e.g. 


Henriksen (20111) 


tand = 


sind' 


7 „(u)(cos0' -I- u/v') ’ 


(HI) 


where u is the relative velocity, v' is the source particle speed, and 7 „ is the bulk flow Lorentz factor. The 
inverse relation is found by interchanging the primes and replacing u by —u. 

This transformation shows that only with u/v' ^ 1 and 7 „ « 1 do the angles coincide and hence that 
isotropy is invariant. Otherwise, even a source particle moving at 9' = 90° when 7 ii = 3 (,5 « 0.94) and 
u/v « 1 is directed at d « 18.4° in the local rest frame. This cone contains all other forward moving 
particles. Particles with 9' > 90° are also beamed to positive angle in the rest frame. 

Consequently, an isotropic distribution of particles in the source frame is a directed beam in the observer 
frame. Fortunately the beaming due to synchrotron radiation itself tends to dominate this effect in the 
integration over solid angle. 

The energy of the particles transforms according to 


E = jM+'^)E/ 


/i u-v', , 

7 = 7«(l H—. 


(B2) 


where v' is the source particle velocity. The inverse follows by interchanging the primes and reversing the 
sign of u. We see that the energy of forwardly moving particles is increased by the roughly constant factor 
2"fu ■ This implies that the power law distribution in energy in the source frame changes only its limits in 
the rest frame. 

It is of interest to clarify the transformation of intensity or ‘brightness’, I,^. this was discussed originally 


in this context in the paper 

Blandford & Kbnigl (1979 

, but it is worth repeating the details 

The radiation pattern of a source transforms as (e.g. 

Henriksen 

(2011)) 


d^E _ I d^E' 
d^dt d^'df' 


(B3) 


where k = (I — k • “) and H denotes solid angle. The unit vector k lies along the direction of a ray. For a 
distant observer however the time is tobs = t + R/c and hence dtobs = dtn, given R as the distance to the 
observer along the emitted ray. Consequently we have 

d'^E _ I d^E' 
diddtobs 7u^"^ dQ,'dt' 

However to calculate the relation between the emissions per unit frequency we must use the Doppler shift 
relation dv = dv'/{’^u) to obtain 


d'^E _ I d'^E' 

dfldtobsdv dQ,'dt'dv' 


(B5) 
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Finally, to obtain the transformation of the intensity, we must consider an area transformation. It is not 
simply given in terms of the Lorentz contraction since the information about the moving surface is delivered by 
the light rays crossing it, rather than from measurements made by two different but simultaneous observers. 
If we imagine a surface element da' oriented perpendicular to the plane of k and u with its normal direction at 
an angle 9' to the relative velocity (hence parallel to a ray in that direction), then the component da' sin O' is 
affected by the change in reference frame. This component becomes equal to da sin 9 in the reference frame 
after tracing the light rays (which define in fact the angle transformation) . Hence we have by equating these 
two expressions that 


da' sin 9 

~T~ = 
da sm 9' 


(B6) 


where we used the standard angle transformation for light rays to get the last expression on the right. In 
the special case of longitudinal motion so that 9' = 9 = 0, this last expression continues to hold (see e.g. 


Ellis fc Williams (2000)). 


We are now in possession of the transformation of intensity in the form 


I ( )= = 1 d^E' ^ 

" dQ.dtobsdvda 7^^^ dVt'dt'dv'da' 7^^^ ’ 


(B7) 


which must be used together with v' = ly^uK- The inverse may be found in the usual way. We see that the 
spectral index is essentially invariant between the two frames over a small range in angle, which is expected 
for synchrotron radiation. The comoving radiation energy density is reduced from that measured in the rest 
frame. The co-moving frequency is reduced by the Doppler factor which might be thought to encourage self¬ 
absorption in the co-moving frame. However the electron number density is also reduced by 7 „. If moreover 
the electric field in the source frame is zero, then the perpendicular magnetic field transformation is given 
by B'j^ = B±/"fu. This reduction also offsets the frequency reduction. In the end the source absorption 
coefficient is virtually invariant (for p = 2) relative to the rest coefficient. 

Another important invariant is the polarization of the source, which we discuss in Sect. 4.4 and Appendix 
C. This invariance is due to the invariance of the phase of the electromagnetic wave. 


In summary, the above relations show that, unless the jet outflow is highly relativistic, i.e. 7 „ >> I, 
measurements in the observer’s rest frame are adequate to describe the instrinsic angles, energy and magnetic 
fields of the source. Note that 7 „ ~ 1 still encompasses the possibility of very rapid jet speeds. For example, 
u = 0.5c => 7 „ = 1.15. 

We have neglected the cosmological red-shift since NGC4845 is so close. However, this affects the observed 
frequency by a well known conversion from the local rest frame frequency (oc 1/(1 -|- z)). 
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C. Spectral Evolution of a Simple Ensemble of Relativistic Electrons 


Let us initially assume that the source is, in part, a stationary ensemble of relativistic electrons that is 
the result of a unique acceleration and that the event that has produced the acceleration occurred at the 
time of the hard X-ray flare at time, to- 

For an individual electron, the synchrotron lifetime, tg, may be written as 


tg = “ 


464 




years, 


(Cl) 


where B±{—2) is the perpendicular (to the electron orbit) magnetic field in units of 10“^ G and ^5 is the peak 
frequency in units of 5 GHz. The critical frequency, essentially the peak frequency produced by electrons of 
a given energy, implies that 

,1/2 


7 e = 345- 


(C2) 


H^(-2)i/2- 

where 7 e is the Lorentz factor of the relativistic electron radiating at 5 GHz. 

For an ensemble of particles in a magnetic field, emission near the spectral turnover will be dominated 
by those electrons of energy, 7 eTOeC^, where rrie is the electron mass and c the speed of light. Thus, from 
Eqn. Cl one expects the peak frequency to decline as (t/to)~^- 


From Table we see that the peak frequencies are in the ratio vt 2 Ivti = 4.03/4.87 = 0.83 after T1 -|- 56 
days, and in the ratio vr’il^Ti = 3.21/4.87 = 0.66 after T1 -|- 196 days. Taking T1 to be 342 days after the 
X-ray peak at to (Sect.|^, the expected ratios are ((Tl -|- 56)/Tl)“^ = 0.74 and ((Tl -|- 196)/T1)“^ = 0.40, 
respectively. The dependence of the peak frequency on time is much closer to being linear than it is to being 
the quadratic variation expected from such a stationary ensemble emitting near the peak of the spectrum. 
For 1^5 = 1 (a spectral turnover of 5 GHz, Fig.|^, for fields even as large as 0.04 G (a comfortable upper 


limit, see Artyukh et al. |[2008 ) the lifetime of particles radiating at the spectral turnover is about 58 years 


and 7 e = 173. Higher energy particles by about a factor of one hundred would be evolving substantially. 
However the corresponding frequency is in the infrared. Allowing for the magnetic field to change over the 
course of the observations is the basis of our explanation of the frequency peak time dependence in the SSA 
model (see Sect. 4.3). However in this case either the or time dependence in magnetic field gives 


quite the wrong time dependence for the peak frequency. Moreover the evolution timescale becomes longer 
with a weakening field which results in an even poorer fit. We observe substantial evolution in less than half 
a year at GHz frequencies, however, which tends to exclude this simple picture even if the peak frequency 
ratios were acceptable. 

We conclude that the behaviour of a source in which there has been a unique acceleration of electrons 
and whose emission is now fading with time does not match the observations. 
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D. Thermal Absorption 


Can thermal absorption explain the low frequency cutoff? The thermal absorption coefficient may be 
written, by ignoring some small logarithmic terms, as (Pacholczyk ||1970| 


n: 


^ 3/2 


( 9)2 


X 10 


-"5 cm-\ 


(Dl) 


where the subscript, ec, refers to ‘cold’ (i.e. non-relativistic, or thermal) electrons, the frequency, v{9), is in 
GHz and the temperature is in units of 10^ K. Multiplying this by the line of sight distance, s, and requiring 
the product to be one (for unity optical depth) shows that, on a scale of 10^^ cm in the optically thick limit 
(^(9) ~ 1), we need a cold plasma density of Uec ~ 10"^ cm~^. This might be in a ‘sheath’ around the source 
(i.e. within or very near the cone). On a scale of one kpc one needs only « 60 cm~^. Such an ‘HII region’ 
might easily intervene along the line of sight in an edge-on galaxy. 

However if we calculate the effective spectral index, Oe//, i.e. the observed spectral index that results 
from a synchrotron spectrum, I^{o) = Ciy°‘, with an additional absorbing thermal screen, then the observed 
spectrum will be 




(D2) 


Differentiating this spectrum with respect to frequency and noting the dependence of (Eqn. Dl) 
results in an observed spectral index of 

tte//= a-|-2 kj^s. (D3) 

For a non-thermal spectral index of a = —0.5 (Table[^, then we calculate for Kec„s = 1 that cteff = 1.5, 
which is too steep (again, Tablej^. At v{9) > 3, the absorption coefficient is down by a factor ten. Therefore 
we obtain Oe// = —0.3, which does not fit the spectra in Fig. Therefore, thermal absorption does not 
describe the turn-over in the spectrum and synchrotron self-absorption (Sect. 4.2.1) is preferred. A more 
detailed analysis based on best fits to the data is given in the text, but our basic conclusion remains. 
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E. Circular and Linear Polarization in the Beckert & Falcke Formulation 


A detailed calculation of the circularly polarized flux from the jet model described in Appendix [A| must 
be left for elsewhere. This calculation has been done in Jones (1988), although frequency dependences 


for the transition region were not given in that paper. However, a quantitative explanation of how linear 
polarization is converted to circular polarization is indeed possible, as will be described here. 

In this section, we adhere closely to the formulation in the appendices of [Beckert fc Falcke (2002). 
We adopt an isotropic relativistic plasma with a power law energy distribution having p = 2 (s in their 
notation), as suggested by our spectral index in the optically thin limit (Table |^. This plasma state implies 
that Ku = Vu = hg = 0 , using the notation of Beckert & Falcke ( 2002[ ) for the U absorption coefflcient, the 
U emissivity, and the extraordinary wave conversion coefficient. Our only additional assumptions are to take 
Kv as negligible and to neglect KqV compared to dPV/ds^. The former assumption is justified by the stonger 
decline with frequency (by a factor \/ vb^_/v) of kv compared to the remaining absorption coefficients. The 
latter assumption neglects the second order change in V by absorption over the whole source, compared to 
a local second order change. 

The ultrarelativistic conversion coefficient given in Sazonov (19691 requires taking a limit as p —> 2 which 
we obtain as 



Kc = 2ko(z/)( —)(sin(^)^ In 

V 

where the scaling opacity (proportional to the Faraday absorption coefficient is 

,,2 


The plasma frequency v„ is 


c V 


2 

^v = - 


2 ’ 


Trme 


(El) 


(E2) 


(E3) 


and 7 o is the minimum Lorentz factor in the distribution. 


A straightforward but tedious manipulation of the equations in Appendix A of Beckert & Falcke (2002) 


using these assumptions, yields the following solution for the Stokes parameters as a function of path length 
in the source, s: 


w = 
= 


KlK^ 

VV 

Qu = —V- 

Kp 


-(kj — Kq + Ki)(l — e H -(k^ —Kj)e ”^'^(1 — cosks + (tan$) sinKs), 

(k^ — /C/)(K/e“'‘^®(sinKs + (tan$) cos ns) — k), 


(E4) 


r]vKi 


{k^ — «:/)(e '^^®(cosKS — (tan$)sinKs) — 1). 


KcKpK^ 

The total intensity may be found from 

nj Kp 


f V, ds+'^iUoe-^^^-U,) + Uoe 

J Kp 


(E5) 


We note that in an optically thick region (k/L large and kq/kp small), « rji/ni, which is proportional 
to as expected (in this appendix, L is the total path length through the source, whereas s is a general 
path length coordinate). 

The initial conditions are H(0) = <5(0) = 0, /(O) = t/(0) = Uo and we have set to zero an arbitrary 
constant of integration that gives a spurious linear increase of V with s. The coefficients, k/, kq = (9/13)k/, 
Kp^ and the emissivity, py, are as in||Beckert fc Falcke (|2002 ). In addition we have defined 


2 / 2 2 I 2 \ 

K — (Kc ~ 


(E6) 
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and 


tan $ = — 

Kj 


KnK 


Uo 


ki{kI - nj) rjv' 


(E7) 


The linear dependence on ryy is a result of isotropy, our boundary conditions and the identity k/? 7 q — 
= 0 j which hold for the coefficients (Beckert & Falcke 2002). However the solution in detail as 


displayed below shows that V is mainly dependent on Uo, as is expected. 

We turn to an analysis of the frequency dependences of these solutions. It is convenient to express 
all frequencies in terms of one GHz (denoted vg), and to let all emission and absorption coefficients have 
their frequency independent values at one GHz when expressed with their usual names () modified to () . 
Thus r]v = fwjvg and Kc = i^i = and so on. Only n requires a slight complication since 

Kp = Kpjv^. Hence from the definition, k may either have the form k = R/vg if Kc, kq are dominant, or 
alternately k = Rpjv^ if Kp is dominant. Other than in these extreme cases, k presents a mixed frequency 
behaviour. 

In the next two subsections, we proceed by assuming that k/L 1 (T is the path length through the 
source) and consider the limits of small and large Faraday depths. 


E.l. The Case of Small Faraday Depth 


In the first case, let us assume that n,p is small compared to or pq. we note that this requires that 
the minimum Lorentz factor 70 be essentially equal to the 7e associated with the peak frequency in the radio 
band. Then the solution for V becomes, to first order in k/L and kL, 


V^tjvL 


ki{kI - «:?) 


(tan$) + r]vL{l + 


')• 


(E8) 


In addition taking np small, the expression for the circularly polarized intensity (cpi) becomes with the 
frequency dependences displayed 


~ r - Kj) tan$ ^ K? - Kq 1 

^ - ~2 \ 3/2 - \-r]vL^ - ^2 -, 

(Kc - Vg Kj . - Kj ^ 1/9 


(E9) 


where 


tan $ = 


^ (Kc - Kq)^/^ 1 


Kl Ijv K/ Rj — Rj Vg 


Gonsequently after inserting this value we have simply 

E 


fjvL UoRcL 




(ElO) 


(Ell) 


Thus when the Faraday absorption is small compared to relativistic conversion, there is an internal 
generation term in the circular polarized flux oc v~^ plus a conversion term cx It is the latter term 

that we are particularly interested in since this is the term that dicates the linear to circular polarization 
and also has the steep spectral behaviour that is observed (Table |^. However, the first term could also 
contribute since the spectral indices, on average, tend to be slightly flatter than -3. 


A numerical calculation using the values given in Beckert & Falcke 
(L(17) is in units of 10^^ cm and 70 = 200) 


(2002) and our formula for Kc gives 


RoL « 7.7 X 10-3L(17)H2(-2)ne(100), 

fjvL « 0.28H2(-2)ne(100)L(17) X 10® {= 22 6*2 H2(-2)ne(100)L(17) mJj/}. (E12) 
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The actual circularly polarized flux measured at the antenna Vi,a = ^sV. Thus the first term in Eqn. Ell 
must be multipl ied b y the factor fig, This factor is « 7.8 x so that the numerical factor in the 

second of Eqns. E12 becomes 22 0^. The source angular size is expressed in milli-arcsec, and may be as 


small as 0.1 when L{17) = 1. For the second term, we will consider Uo as if it has already been corrected to 
flux density, thus Eqns. [ETT and E12| together give 


V 


22 0/ B^(-2)ne(100)L(17) 7.7 x 10-^UoL{17)B^{-2)neil00) 




+ 


mJy. 


(E13) 


Hence the first term in V may contribute to our circular polarization if B{—2) « 4 and the relativistic 
electron density is as large as 100cm“^ at 1.5 GHz {V = 2.3 mJy). Under these conditions the 

second term would require a linearly polarized flux Uo ~ 60 mJy, in order to produce the observed circularly 
polarized flux by conversion (V = 2.2 mJy), i.e. comparable values that are roughly consistent with our 
total circularly polarized flux density. It is the second term that is favoured by our measured in-band values 
of ay however. The required incident Uo that makes the second term comparable to the first term is 
larger than any linearly polarized flux (close to 24 %) that we or others observe. We must then consider 
the exit values of U and Q further below, especially after passing through a possible depolarizing screen (cf. 
O’Sullivan et al. [2013 1. 

The dominant term in the total intensity is approximately the first term in Eqn. |E5[ Assuming that both 
terms in Eqn. |ElT] are comparable, we obtain a rough estimate of the percentage cpf as 


Vv 




— cot ( 


I b^{-2) 

vg 


X 0.9 X 10“ 


(E14) 


This can be of order 2 % for B{—2) > 4. 

We note that any balanced mixture of the two terms in Eqn. |E11| will change quickly with increasing 
frequency towards a behaviour. Moreover the cpi can not be much larger than KcLUo (cf. O’Sullivan 
et al. |[2013 ). Because of the harmonic functions of ks in the general expression for V, this cpi state also 
changes rapidly in the optical depth transition zone as we require. 

Unfortunately, proceeding in the same fashion as for the circular polarization, we find that the Stokes U 
is modified only by absorption. It becomes 


U = Uo{l-'^), 


(E15) 


and so is essentially unchanged from Uo. The latter quantity is thus to be constrained by the small value of 
linearly polarized flux (Ipf) that we observe. One must appeal to a depolarizing screen to remove the linear 
polarization (e.g. O’Sullivan et al. 1|2013 | if it is really of this order. We note that a magnetic field as large 
as B{—2) = 10 (e.g. Artyukh et al. ||2008 ) would reduce the required Uo to a more reasonable 3 — 4%, but 
this would still have to be removed by an inhomogeneous, depolarizing screen to satisfy our limit. 

Finally we calculate Stokes Q in the same approximation and limit. When there is very little Faraday 
rotation we find that Q « 0 in lowest order provided that k » Kc , and that this quantity is much larger 
than either ki or Rq. 


E.2. The Case of Large Faraday Depth 

A second limiting case arises when the Faraday depth is large so that Kp is large compared to Kc or kq 
(and therefore k/ since for p = 2 this coefficient is larger only by the factor 13/9) . In this limit k ■ 
and we find at small kjL that 




V : 


VyL 

vg 


_ kUk'^-k'j) tan4) , , 

W -— -4 - - - ^ sm (kfL). 


(E16) 
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Moreover 


^ Kf Kc Kl 

tan $ = Fg + — — 

Kj Ki — Kj r]v 


Uo 

1 ~ ^9 ’ 


and hence 




riv hi - hi 


sin {npL) + 


fivL Uohc 


hpFg 


sin (KpL). 


(E17) 


(E18) 


l^g rvp Vg 

The terms in fjy must be assumed to be reduced numerically by Us to give the observed flux, as in the 
previous case. 

The first term in Eqn. |E18| is small compared to the second term around Fg = 1 according to our limiting 
assumption. Thus the combined frequency dependence due to the second and third terms is not changed 
from the small Faraday rotation limit. The magnitude of the conversion term (proportional to Uo) is however 
reduced by the factor hcjhp, which is assumed to be quite small. In fact this ratio is (see formulae in Beckert 
k Falcke |[2002l ) 


hr 

— « 0.75 X 

np V 


sin 




)■ 


, X- . y (E19) 

cos (j) 2 In 7o 

This is very small unless 7o approaches 100 or more, which becomes the previous case. Hence in this case 


the conversion term is not important compared to the synchrotron emission term given in Eqn. E12 This 
can give sufficient cpf as we have seen, and it is not oscillatory. However the frequency dependence is 
rather than the behaviour we observe. 

For completeness, and in order to add insight into the effects of Faraday inhomogeneity, we calculate the 
Stokes parameters Q and U in the Faraday dominant limit as 


and 


U FiUo cos (kpL), 


Q = —UoSui{KpL). 


(E20) 


(E21) 


The only way to suppress a large Uo is to imagine many cells with random properties. Then L is replaced 
by a characteristic s for each cell and, in the mean, the linear polarization is suppressed Beckert k Falcke 


(2002). The same suppression would occur in a depolarizing sheath. 
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F. X-ray Emission via Inverse Compton 

Appendix [X| outlines the jet model in which we imagine that the X-ray emission around the time of the 
peak of the X-ray light curve, comes from a small X-ray ‘lobe’ (see Fig. [^. We wish to link the currently 
observed radio emission to a time approximately a year earlier to determine whether conditions at that 
time could account for the X-ray peak flux via inverse Compton upscattering, also known as Synchrotron 
Self-Compton, or SSC. As before, our calculations are order-of-magnitude to check for feasibility. 

As indicated in the text (Eqn. C2), the critical frequency, Vc (Hz), for a relativistic electron of Lorentz 
factor, 7e, in a magnetic held of strength, B, is 


= 4.2 X 10'^7e^Hj_(-2) 


(FI) 


where B±{—2) is the perpendicular magnetic held strength {B± = Bsind, with 9 the pitch angle). Hence 
for an electron to radiate in the radio regime at 5 GHz (the peak of the spectrum at time Tl, Table [^, 
in a magnetic held of B±{—2) = 1, a Lorentz factor of 7e « 347 is required. Such a value could describe 
conditions in the radio source around the time of the radio observations. 

We need to consider the earlier time in the X-ray lobe at which 7e and B±{—2) are both higher and the 
source size is smaller. The SSC requirement is i^x = 7e^ so from the above equation. 


ux — 4.2 X 


10Se^Hj_(-2) 


(F2) 


The hard X-ray light curve of Nikolajuk & Walter (2013) applies to the energy range, 17.3 - 80 keV 
and we adopt 30 keV, or vx = 7.2 x 10^*^ Hz as the X-ray frequency. Also, for a jet of constant velocity, 
Bj. (X (see Appendix®, so for an X-ray lobe that is a factor of 10 smaller than the current radio source 
size (10^® rather than 10~cm, or 0.1 mas rather than 1 mas), the magnetic field strength can be estimated 
to be i?x(— 2) = 100 at the earlier time. We can now find the Lorentz factor for relativistic electrons in 
the X-ray lobe, which from Eqn. F2 is 7e = 1109. Then, by Eqn. FI the synchrotron critical frequency 
would be Vc = 5.5 x 10^^ Hz or Ac = 55 /rm. Thus the synchrotron emission in the X-ray lobe would be 
in the infra-red at the earlier time. Though not required for this development, one can show via standard 
formulae, that such emission would be optically thin. 

We now need to calculate the infra-red flux in the X-ray lobe, given these conditions. First, we find the 
specific intensity of the synchrotron emission, specifying infra-red with the subscript, IR, 


= e. 


sx 


(F3) 




where is the emissivity and sx is the line-of-sight distance through the X-ray lobe. Here, C^{p) and 


Cl are Pacholczyk’s constants (Pacholczyk 1970) for electron energy spectral index, p, and is a 


normalizing factor for the energy spectrum, (cf. Sect. 4.3). We have explicitly expressed the magnetic field. 


Bjr separated from sinOm because, in the next step, we will compare current radio values to IR values at 
the earlier time and, when 7e changes, so does the pitch angle. 


We now multiply and divide Eqn. F3 by the same expression for the (current) radio emission, Ivr, where 


R refers to the radio, finding after some manipulation. 




= IVB. 


OIR 


OR 


Bir \ 

Br J 
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/ sin{9iR) \ 
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(l-p)/2 


■Sx 
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(F4) 


A conversion to flux. Sir, requires multiplying by the source solid angle and by the frequency. If we again 
take p = 2 which is implied from the frequency spectral index in the radio, then Eqn. |F4| becomes 


3/2 


3/2 


- 1/2 


c c , Noir\ (BirY^ fsin{9iR)Y'^ f^iRV ' fsx\ 


(F5) 
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We can now relate the radio quantities, which are mostly known, to the infra-red quantities at the earlier 
time to estimate Sm. We take Sr = vr where vr = h x 10® Hz and = 433 mJy (Table [ t] fo r time 
Tl) to find Sr = 2.1 x 10“^"^ erg s“^ cm“®. For an adiabatically expanding source, we found (Sect. 4.31, that 
Nqjr/Nqr oc where tR, tjR, are the times for the radio emission and IR emission, respectively. 

We can take Iir/Ir « 1/2, since the time from the outburst to the radio observations is about twice the time 
from the outburst to the X-ray peak, leading to Noir/Nqr « 1.6. The ratio, Bjr/Br = 100, from above. 
For the pitch angle, we take sin{6j r)/ sin{9 r) = "ieR/leiR = 347/1109 = 0.31. The ratio of frequencies is 
virIvr = (5.5 X 10^®)/(5 X 10®) = 1.1 X 10® and the ratio of sizes is sxIsr = 0.1. 

The result is Sjr = 1.9 x 10“®® erg s“® cm“^. 


It remains to compute the X-ray flux from SSC. From Eqns. 23 and 24 expressed in a more convenient 
form in terms of the infra-red synchrotron emission in the shocked region that emits the X-rays, 


^Xssc 


= s 


IR 


Bradi^ 

u, 


= 5.7 X 10 


12 




mag 


9,^Bir{-2Y 


(F6) 


for Sxssc’ ^iR s“® cm“®, 9s in mas and Bir{—2) in units of 10“^ G. With the above value of Sir, 

9s = 0.1 and Bir{—2) = 100 (see above), we find, Uradju/Umag = 1-1 and Sxssc — 2-0 x 10“®® erg s“® 
cm“^. The peak of the X-ray light curve is Sx ~ 8.0 x 10“®® erg s“® cm“®. 

Thus, our estimated X-ray flux from SSC falls short of the observed X-ray peak by only a factor of 4. The 
result is quite sensitive to the choice of magnetic field in the radio which was taken as Br{—2) = 1 above. 
For example, if Br{—2) = 4 then Uradm/Umag = 0.06 and Sxssc — 1-0 ^ 10“®^ erg s“® cm“® which falls 
short of explaining the X-ray peak. If Br{—2) = 0.5, then Sxgsc — 1-0 ^ 10“®® erg s“® cm“® which is 
more than adequate to explain the X-ray peak; however, Uradin/Umag = 4.6 which is uncomfortably large. 

Since the result is so sensitive to the magnetic field, one needs to examine it more closely. We have 
assumed that B^ oc r“^ which accounts only for geometric expansion. Yet, our jet model indicates that the 
X-ray lobe is a shocked region that should enhance the magnetic field. If we allow for geometric expansion, 
as before, but include a modest enhancement to the magnetic field in the X-ray lobe of a factor of 1.8 (a 
strong shock would give factor of 4) then Uradia/Umag = 1-1 and Sxssc = 8 x 10“®® erg s“® cm“®. With 
these conditions, SSC can account for the X-ray flux at the peak of the light curve. 


Facilities: VLA. 
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